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ABSTRACT. In 1979, building on S. Lie's theory of symmetries of (partial) differential equa- 
tions, RJ. Olver formulated inductive formulas which are appropriate for the computation of the 
prolongations of an infinitesimal Lie symmetry to jet spaces, for an arbitrary number n > 1 
of independent variables (x 1 , . . . ,x n ) and for an arbitrary number m > 1 of dependent vari- 
ables (y 1 , . . . , y m ). This paper is devoted to elaborate a formalism based on multiple Kronecker 
symbols which enables one to handle these "unmanageable" prolongations and to discover the 
underlying complicated combinatorics. Proceeding progressively, we write down closed explicit 
formulas in four cases: n = m = 1; n > 1, m = 1; n = 1, m > 1; general case n > 1, 
m > 1. As a subpart of the obtained formulas, we recover four possible versions of the (multi- 
variate) Faa di Bruno formula. We do not employ the classical formalism based on the symmetric 
algebra (cf. e.g. H. Federer's book, p. 222), because it hides several explicit sums in symbolic com- 
pactifications, and because the presence of supplementary complexities (e.g. splitting of indices, 
combinatorics of partial derivatives) impedes us to apply such compactifications coherently. Our 
method of exposition is inductive: we conduct our reasonings by analyzing several thoroughly or- 
ganized formulas, by comparing them together and by "drifting" towards generality, in homology 
with the classical style of L. Euler. 
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§1. Jet spaces and prolongations 

1.1. Choice of notations for the jet space variables. Let K = K or C. Let n > 1 and m > 1 

be two positive integers and consider two sets of variables x = [x , . . . , x n ) £ IK™ and y — 
(y 1 , . . . , y m ). In the classical theory of Lie symmetries of partial differential equations, one 
considers certain differential systems whose (local) solutions should be mappings of the form 
y = y(x). We refer to IOL1986I and to |BK1989| for an exposition of the fundamentals of 
the theory. Accordingly, the variables x are usually called independent, whereas the variables y 
are called dependent. Not to enter in subtle regularity considerations (as in |M2005|), we shall 
assume C°°-smoothness of all functions throughout this paper. 

Let k > 1 be a positive integer. For us, in a very concrete way (without fiber bundles), the 
HL-thjet space J^ m consists of the space K" +m+m ^. m \ equipped with the affine coordinates 

d- 2 ) ( x \y 3 MiM^ >»ii,fa,.--,*« 

having the symmetries 
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for every A with 1 < A < k and for every permutation a of the set {1, . . . , A}. The 
variable y? { i is an independent coordinate corresponding to the A-th partial derivative 

dx i 1 Q x i^ dx i x , which explains the symmetries (1.3). 
X In the classical Lie theory ( IIOL 19791 . IOL 19861 . IBK1989I ). all the geometric objects: point 
transformations, vector fields, etc., are local, defined in a neighborhood of some point lying in 
some affine space K N . However, in this paper, the original geometric motivations are rapidly 
forgotten in order to focus on combinatorial considerations. Thus, to simplify the presentation, 
we shall not introduce any special notation to speak of certain local open subsets of K n+m , or of 
the jet space J£ m = K n+m+m m m < ( etc.: we will always work in global affine spaces K N . 

1.4. Prolongation ip^ of a local diffeomorphism ip to the n-th jet space. In this paragraph, 
we recall how the prolongation of a diffeomorphism to the K-th jet space is defined ( IOL 19791 . 
IQTT3561 . IBK1989I). 

Let G K™ be a central fixed point and let ip : K n+m — > K n+m be a diffeomorphism whose 
Jacobian matrix is close to the identity matrix at least in a small neighborhood of x„ . Let 

(1-5) J x , := i K x *^vii 1 ^vii 1 ,i 2 ^ 'y*ii,ig,...,inj ^ ^n,m\ Xt 

be an arbitrary ft -jet based at x*. The goal is to defined its transformation tpW (J x ) by tp. 

To this aim, choose an arbitrary mapping K™ 3 x i— > g(x) G K m defined at least in a 
neighborhood of x* and representing this K-jet, i.e. satisfying 

d-6) d u ....M - d J. g3 . dx iS x *)' 

for every A G N with < A < k, for all indices with 1 < i\, . . . , i\ < n and for every 

j £ N with 1 < j < m. In accordance with the splitting (x, y) G K" x K m of coordinates, 
split the components of the diffeomorphism cp as cp = (4>,tp) G K™ x K m . Write (x,y) the 
coordinates in the target space, so that the diffeomorphism ip is: 

(1.7) K™+™ 9 (x,y) ^ (x,y) = (<t>(x,y),^(x,y)) G K n+m . 

Restrict the variables (x,y) to belong to the graph of g, namely put y := g(x) above, which 
yields 



(1-8) 



x = (f)(x,g(x)), 
V = ^(x,g{x)). 



As the differential of (p at is close to the identity, the first family of n scalar equations may 
be solved with respect to x, by means of the implicit function theorem. Denote x = x(x) the 
resulting mapping, satisfying by definition 

(1.9) x = cp{x{x),g(x{x)))- 

Replace x by x(^) m tne second family of m scalar equations (1.8) above, which yields: 

(i-io) y = i>(x(x),g(x(m- 

Denote simply by y = ~g(x) this last relation, where ~g(-) := ip (x(')i 5(x( - )))- 

In summary, the graph y = g(x) has been transformed to the graph y = ~g(x) by the diffeo- 
morphism ip whose first order approximation is close to the identity. 

Define then the transformed jet ip( K > (J£,) t0 b e me K-fhjet ofg at the point x* := </>(x*), 
namely: 

(lid ^( J xJ ; = iFz ^n e J Z*\* ■ 

\OX ■■■OX / !<»! i\<n, 0<X<k 
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It may be shown that this jet does not depend on the choice of a local graph y = g(x) representing 
the K-jet at x*. Furthermore, if ir K := J* m — ► IK" 1 denotes the canonical projection onto 
the first factor, the following diagram commutes: 



7 K 

<-'n,T 



rn-\-m 



1.12. Inductive formulas for the K-th prolongation tp^ K K To present them, we change our 
notations. Instead of (x,y), as coordinates in the target space K' 1 x K m , we shall use capital 
letters: 

(1.13) (X\...,X n ,Y\...,Y m ). 

In the source space K n+m equipped with the coordinates (x, y), we use the jet coordinates (1.2) 
on the associated K-th jet space. In the target space K n+m equipped with the coordinates (X, Y), 
we use the coordinates 



(1.14) 



(x i V" V Y j Y j ) 



on the associated K-th jet space. In these notations, the diffeomorphism tp whose first order 
approximation is close to the identity mapping in a neighborhood of x % may be written under the 
form: 



(1.15) 



tp : (z*V) h- (X\Y 3 ) = (XW ,yt'), K 3 '(^'./)) . 



for some C°° -smooth functions X % (x l , y 3 '), i — 1, . . . , n, and Y- i {x l ,y° ), j = 1, . 
first prolongation ip^ of tp may be written under the form: 



, m. The 



(1.16) 



x l , y 3 , yl 



X\x\y 3 ), YHx\yi), Y xH [x l ,y> ,y>, 
which depend on the pure first jet variables y J , . The way 



for some functions YL 

how these functions depend on the first order partial derivatives functions X" 1 ,, , X 1 ., , Y 3 ., , Y 3 ., 

r ± x z yj 1 x i yJ 

and on the pure first jet variables y 3 ., is provided (in principle) by the following compact formulas 
( IBK1989I ): 



(1-17) 



V yL 



DlX 1 ■ 


■ D\X n \ 




( d\y ] \ 




■ DiX n ) 







where, for i' — 1, . . . , n, the Dj, denote the i'-th first order total differentiation operators: 
(1.18) 



d ill r\ 
\ - ,•' O 



J'=l 



dyi 



Striclty speaking, these formulas (1.17) are not explicit, because an inverse matrix is in- 
volved and because the terms D\,X l , D\,Y 3 are not developed. However, it would be ele- 
mentary to write down the corresponding totally explicit complete formulas for the functions 

. Y i : (•'"•// / •///•'). 

Next, the second prolongation tp^ is of the form 
(1-19) 



(>> (V,/4) , Y xiixi2 (^,^',4,4^)) . 
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for some functions Y xHxi2 (x l , y J , yj, , y?, ., ) which depend on the pure first and second jet 



l 1 &l,i 2 



variables. For i = 1, . . . ,n, the expressions of Y xiixi are given by the following compact 
formulas (again IBK1989I ): 



(1.20) 





IX 1 1 


-( 









DlX 1 



D\X' 



DlX" 



n 2 Y° 



D 2 Y 3 
^n 1 xn 



where, for i' = 1, . . . , n, the D 2 denote the i'-th second order total differentiation operators: 



(1.21) 



Dl 



d 

dx 1 ' 



Ed 



d 



m n 



y y J ., _iL 



j'=i y j'=i i{=i 



Again, these formulas (1.20) are not explicit in the sense that an inverse matrix is involved and 
the terms D\,X l , D 2 Y xil are not developed. It would already be a nontrivial computational task 
to develope these expressions and to find out nice satisfying combinatorial formulas. 

In order to present the general inductive non-explicit formulas for the computation of the n- 
th prolongation tp( K >, we need some more notation. Let A £ N be an arbitrary integer. For 
i' = 1, . . . , n, let D\, denotes the i'-th X-th order total differentiation operators, defined precisely 
by: 
(1.22) 



D 



A .. 



d 
dx 1 ' 







m n 



1 X X fop + X X y i'AA ^ 



m n 



j'=l ti,»i,...,»i_l=l 



j' = l i[,i' 2 = l 

d 



1 ,l 2 ,...,l x l 

i' ,%',,...,% 



dy i , 



Then, for i = 1, . . . , n, the expressions of Y^. 4l x i A _ i: ^ 
formulas (again IBK1989I ): 



/ Y j . 

x % i—x'^- l x 



(1.23) 



Y 3 . 

Xn...x^-iJf 



DlX 1 



V ^ • • ■ D n XU J 

Again, these inductive formulas are incomplete and unsatisfactory. 

Problem 1.24. Find totally explicit complete formulas for the n-th prolongation 



are given by the following compact 

d\x- y 1 ( d^y xii x%x _ 1 ^ 



DiYi 



X*i...X'A-i 



Except in the cases n = 1 , 2, we have not been able to solve this problem. The case k = 1 is 
elementary. Complete formulas in the particular cases k = 2, n = 1, m > 1 and n > 1, m = 1 
are implicitely provided in I M2004"a) and in |M2004bj|, where one observes the appearance of 
some modifications of the Jacobian determinant of the diffeomorphism (p, inserted in a clearly 
understandable combinatorics. In fact, there is a nice dictionary between the formulas for ip^ 
and the formulas for the second prolongation C^ 1 of a vector field C which were written in equa- 
tion (43) of 1GM2003I {see also equations (2.6), (3.20), (4.6) and (5.3) in the next paragraphs). 
In the passage from (p^ to C^ 2 \ a sort of formal first order linearization may be observed and 
the reverse passage may be easily guessed. However, for k > 3, the formulas for ip^ K > explode 
faster than the formulas for the K-th prolongation of a vector field C. Also, the dictionary 
between tpW and disappears. In fact, to elaborate an appropriate dictionary, we believe that 
one should introduce before a sort of formal (k — 1) order linearizations of <p( K \ finer than the 
first order linearization CS K >. To be optimistic, we believe that the final answer to Problem 1.24 
is accessible. 
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The present article is devoted to present totally explicit complete formulas for the K-th pro- 
longation of a vector field £ to J* m , for n > 1 arbitrary, for m > 1 arbitrary and for k, > 1 
arbitrary. 



1.25. Prolongation of a vector field to the n-th jet space. Consider a vector field 
(1.26) 



do; 4 

1=1 3=1 



defined in K n+m . Its flow: 

(1.27) (f t (x, y) := exp (tC) (x, y) 

constitutes a one-parameter of diffeomorphisms of K n+m close to the identity. The lift {^tp K ' to 
the K-th jet space constitutes a one-parameter family of diffeomorphisms of J^m- By definition, 
the n-th prolongation D- K > of £ to the jet space J£ m is the infinitesimal generator of (ipt)^ K \ 
namely: 

d 



(1.28) 



t=o 



(<Pt) 



1.29. Inductive formulas for the K-th prolongation As a vector field defined in 

K" +m+m «i mi' , the K-th prolongation may be written under the general form: 



(1.30) 



i=i j 

rn n 



^ ^ 11 dv 3 ^ ^ 11,12 dv - 



+ 



3=1 ti=l 
m n 



E E n, 

j=l ti,...,i„=l 



Here, the coefficients Y? , Y? ■ 



, 4 j are uniquely determined in terms of par- 
tial derivatives of the coefficients X % and y of the original vector field C, together with the 



pure jet variabl es [y 3 ^, ■ ■ ■ , yj lt ___ j K j > by means of the following fundamental inductive formu- 
las ( IOL1979I . IOL1986I . IBK1989I 1: 



(1.31) 



Yl:=^ (y)-J2Di (X k ) yi, 

n 

<<, ■■ (n)-E^ (**) vL,» 



k=l 



Y J - ■ 

* U,»2,— J*<i 



£> K (Y J 

1 tl,»2) — l*K-I 



E (* fe ) ^ 



j 2 ,...,'t/ t -l,fc' 



fe=l 



where, for every A G N with < A < n, and for every i G N with 1 < i < n, the z-th A-th order 
total differentiation operator was defined in (1.22) above. 

Problem 1.32. Applying these inductive formulas, find totally explicit complete formulas for the 
n-th prolongation 



The present article is devoted to provide all the desired formulas. 
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1.33. Inductive methodology. We have the intention of presenting our results in a purely induc- 
tive style, based on several thorough visual comparisons between massive formulas which will 
be written and commented in four different cases: 

(i) n = 1 and m = 1; n > 1 arbitrary; 

(ii) n > 1 and m = 1; n > 1 arbitrary; 

(iii) n = 1 and m > 1; k > 1 arbitrary; 

(iv) general case: n > 1 and m > 1; K > 1 arbitrary. 

Accordingly, we shall particularize and slightly lighten our notations in each of the three 
(preliminary) cases (i) [Section 2], (ii) [Section 3] and (iii) [Section 4]. 



§2. One independent variable and one dependent variable 

2.1. Simplified adapted notations. Assume n = 1 and m = 1, let n S N with n > 1 and 
simply denote the jet variables by: 

(2-2) (x,y,y 1 ,y 2 , ■ ■ ■ ,y K ) £ J£ v 

The K-th prolongation of a vector field will be denoted by: 



(2.3) 



d_ 

dx 



ay oyi oy2 



d 
dy K 



The coefficients Yi, Y2, . 
(2.4) 

where, for 1 < A < k: 
(2.5) D x 



Y K are computed by means of the inductive formulas: 

f Yi := D 1 (y) - D l (X) 2/1, 
Y 2 :=D 2 (Y 1 )-D 1 (X)y 2 , 



Y K :=D K (Y K _ 1 )-D 1 (X)y K , 



dx ^ 1 dy 



V2 



9yi 



dyx- 



By direct elementary computations, for n = 1 and for k = 2, we obtain the following two very 
classical formulas : 



(2.6) 



Yx = y x + [y y - x x \ y\ + [-y y ] ( yi ) 2 , 

y 2 = y X 2 + [2y xy - x x2 ] yi + [y y , - 2 x xy ] { yi f + [-x y2 ] ( yi ) 3 - 

+ [y y -2X x ] y 2 + [-3X y ] yi y 2 . 



Our main objective is to devise the general combinatorics. Thus, to attain this aim, we have to 
achieve patiently formal computations of the next coefficients Y3, Y4 and Y5. We systemat- 
ically use parentheses [•] to single out every coefficient of the polynomials Y3, Y4 and Y5 in 
the pure jet variables yi, y 2 , 2/3, yA and j/5, putting every sign inside these parentheses. We al- 
ways put the monomials in the pure jet variables yi,y 2 ,y3, y4 and y$ after the parentheses. For 
completeness, let us provide the intermediate computation of the third coefficient Y3. In detail: 

Y 3 = D 3 (Y 2 )-D 1 (X)y 3 

+ [y y 2 — 2 X xy \ ( yi ) 2 + [-X yi ] { yi f + [y y -2X X ] y 2 + [-3X V ] y\ y 2 ) 
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y^rj-j + [iy x 2 v - X x s] + [j^y 2 - 2X x 2 y ] (^i) 2 rjj + [-^(tj) 2 ] (yi) j-j-j- 



(2.7) 



] jfej-gj + [-3 x xy ] ymp ^ + [34? 







f [2y xy 2 - X x2y ] (yi) 2 ^+ [y v 3 -2X xy ,] (yi) 3 ^+ [-X y3 ] (yi) 4 ^- 

+ [3V -2^] yijte ryj + [-3^2] ( yi ) 2 y 2 ^ + [2y xy - ^2] j/2 |Yj+ 

+ [jy - 2 2 + [-■*>] 3(yi) 2 y 2 ^ + [-3 *„] (2/2)^ + 



+ [y y - 2 A'x] 7/31-^-1 + [-3 ^] t/it/3 



- [X x ] - [Xy] yn 







(2.8) 



We have underlined all the terms with a number appended. Each number refers to the order of 
appearance of the terms in the final simplified expression of Y3, also written in | BK1989 1 with 
different notations: 

' y 3 = y x s + [3y X 2 y - x xS ] yi + [3y xy 2 -3x X 2 y ) ( m ) 2 + 

+ [y v s -sx xy2 ] ( yi ) 3 + [-x y3 ] ( yi ) 4 + [sy xy -3x x2 ] y2 + 

+ [3y y 2 -9X xy ] Vl y 2 + [-6X y2 ] ( Vl ) 2 y2 + l-3X y ] (y 2 ) 2 + 
+ [y y -3X x } y3 + [-4Xy} ym . 

After similar manual computations, the intermediate details of which we will not copy in this 
Latex file, we get the desired expressions of Y4 and of Y5. Firstly: 

' y 4 = y X 4 + [4y x , y - x xi ] yi + [ey X 2 y 2 -4X x3y ] ( yi f+ 

+ [4y xy3 -6X X 2 V 2] ( yi f + [y yi -4x xy3 ] ( yi ) 4 + [-x yi ] ( yi f+ 

+ [6y x z y -4X x3 ] j/2 + [12^2 - 18X x2y ] yim + 
+ [6^3-24^2] ( yi ) 2 y 2 + [-10X y3 ] (yify 2 + 
+ [3^2 - 12X xy ] (y 2 f + [-15 Xj] y i(y 2 f+ 
+ [4y xy -6X x 2]y 3 + [4y y 2 - 16X xy ) 2/12/3 + [-WX y2 ] (yi) 2 y 3 + 
+ [-10 Xy] y 2 y3 + [y y - 4 X x ] y 4 + [-5 X y ] 2/12/4. 



(2.9) 



Secondly: 



(2.10) 



y 5 = y x s + [5y x * y - x x5 ] yi + [ioy x3y 2 ~5x xiy ] ( y 



+ [10^2^3-10^2 



(mY +[5y xy +-WX x 2 y3 ] (yi) 4 + 

+ [y y5 -5X xyi ) ( yi f + [-x y5 ] ( yi f + [ioy x3y -5x xi ] y2 + 
+ [soy x 2 y 2 - 30X xSy ] yiy2 + [3oy xy3 - mx X 2 y2 ] ( V i) 2 y2 + 

+ [Wy yi -50X xy3 ] ( m f y2 + [-15X yi ] (yi) 4 y 2 + 

+ [15 3V - 30X„, V ] ( V2 ) 2 + [I5y y3 - 75X xy2 ) V i( y2 ) 2 + 

+ [-45X y3 ] (yi) 2 ( y2 ) 2 + [-15 Xy,] ( y2 f+ 

+ [wy X 2 y - wx x3 ) y3 + [2oy xy2 - 4ox X 2 y ) yiy3 + 
+ [10^,3-50^2] ( yi ) 2 m + [-20X y3 ] ( y if m + 

+ [I0y y 2 - 50 X xy ] j/22/3 + [-60 X y2 ] 2/12/22/3 + [-10 Xy] (y 3 ) 2 + 
+ [5y xy - 10 X x2 ] 2/4+ [5 3V -25X xy ] 2/12/4+ [-15^2] (yi) 2 ?/4+ 

+ [-15 Xy] 2/22/4 + [y V -5 Xy] V5 + [-6 Xy] J/ 5 . 
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2.11. Formal inspection, formal intuition and formal induction. Now, we have to comment 
these formulas. We have written in length the five polynomials Yi, Y2, Y3, Y4 and Y5 in the 
pure jet variables 2/1, 2/2,2/3, 2/4 and 2/5- Except the first "constant" term y x *, all the monomials 
in the expression of Y K are of the general form 

(2.12) (yx i r(yx 2 r---(yx d r d , 

for some positive integer d > 1, for some collection of strictly increasing jets indices: 

(2.13) 1 < Ai < A 2 < • • • < A d < k, 

and for some positive integers fii, . . . , /x<i > 1. This and the next combinatorial facts may be 
confirmed by reading the formulas giving Yi, Y2, Y3, Y4 and Y5. It follows that the integer d 
satisfies the inequality d < K + 1. To include the first "constant" term 34«, we shall make the 
convention that putting d = in the monomial (2.12) yields the constant term 1. 

Furthermore, by inspecting the formulas giving Yi, Y 2 , Y3, Y4 and Y5, we see that the 
following inequality should be satisfied: 



(2.14) 



/tiAi + /i 2 A 2 H h MdAd < k + 1. 



For instance, in the expression of Y4, the two monomials (j/i) 3 y 2 and y\ (y 2 ) 2 do appear, but the 
two monomials (2/i) 4 2/2 and (yi) 2 (2/2) 2 cannot appear. All coefficients of the pure jet monomials 
are of the general form: 



(2.15) 



[Ay x o 



B Xy, 



for some nonnegative integers A,B,a,f3 £ N. Sometimes A is zero, but B is zero only for 
the (constant, with respect to pure jet variables) term 34 K - Importantly, X is differentiated once 
more with respect to x and y is differentiated once more with respect to y. Again, this may be 
confirmed by reading all the terms in the formulas for Yi, Y 2 , Y3, Y4 and Y5. 

In addition, we claim that there is a link between the couple (a, (3) and the collection 
{fi%, Ai, . . . , /j,d, Ad}. To discover it, let us write some of the monomials appearing in the ex- 
pressions of Y4 (first column) and of Y5 (second column), for instance: 



(2.16) 



[6y x 2 y 2 
\i2y x 



4x x3y \ ( yi y, 

r xy * ■ 18X x 2 y ] J/12/2, 

[-10 A;.] {yif V 2, 

[±y y * - lQX X y\ J/lJ/3, 



[5 3>x» 



10 P£ X 2y3 



( y 1 



-10 X y ,] ( yi ) 2 y 3 , 



[30^ a3 -60^ v ] ( yi ) 2 y 2 , 
[-15^4] { Vl ) A y 2 , 
[10^-50^] y 2 y 3 , 
[-60^2] j/ij/22/3- 



After some reflection, we discover the hidden intuitive rule: the partial derivatives of y and of X 
associated with the monomial (i/aJ^ 1 • ■ ■ {y\ d Y d are, respectively: 



(2.17) 



x„ 



-Ml^l" 
-H 1 X 1 - 



-^d X d y»l+"-+fd 1 
-^d x d + x J/flH hf-'d 



This may be checked on each of the 10 examples (2.16) above. 

Now that we have explored and discovered the combinatorics of the pure jet monomials, of the 
partial derivatives and of the complete sum giving Y K , we may express that it is of the following 
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general form: 



(2.18) 



K + l 



*-+E E E E 

d— 1 1< Ai <■ ■ ■ < Arf <k /^i >1,...,/^^>1 /^iAiH hMdAd<^+l 

/i(/ti,Ai),...,(/id,Xi) -n _ 

J / x k- i i 1 X 1 Cd*d 1 iM+-+l'd 



n(Mi,Ai),...,(Md,A d ) 



f'd A ci + 1 ^Mei -1 

(yA 1 r---( 2 /A d r. 



Here, we separate the first term y x * from the general sum; it is the constant term in Y„, 
which is a polynomial with respect to the jet variables y\. In this general formula, the 

only remaining unknowns are the nonnegative integer coefficients Ait 1 ' ^ fJ ' d ' g anc [ 
B (/*i,Ai),...,(Md,A«i) £ R In Section 3 beloWj we shall eX pi ain how 

we have discovered their exact 

value. 

At present, even if we are unable to devise their explicit expression, we may observe that the 
value of the special integer coefficients A^ 1 ' 1 ^ and B^ 1 ' 1 * 1 which are attached to the monomials 
ct., j/i, (yi) 2 , (yi) 3 , (yi) 4 and (yi) 5 are simple. Indeed, by inspecting the first terms in the 
expressions of Yi, Y2, Y3, Y4 and Y5, we of course recognize the binomial coefficients. In 
general: 

Lemma 2.19. For k > 1. 



(2.20) 



A=l 

+ [— X y n\ (yi) K + remainder, 



A- 1 



(yi) } 



where the term remainder collects all remaining monomials in the pure jet variables. 
In addition, let us remind what we have observed and used in a previous co-signed work. 



Lemma 2.21. ( IGM2003I . p. 536) For n > 4, nine among the monomials of Y K are of the 

following general form: 



(2.22) 



f Y K = y x n + [Cl y x *-ly - X x n] Vl + [Cl y xK -2 y - C\ X x «-l] y2 - 

+ [cl y x *y - cl x x3 ] Vk _ 2 + [cl y xy - cl x x ,] Vk _ 1+ 

+ [ G l y v 2 - k 2 x xy] ViVk-i + [-Cl X y ] y 2 y K -i+ 
+ [y y - C\ X x ] + [-C\ +l X y ] y x y K + remainder, 



where the term remainder denotes all the remaining monomials, and where C^ := ^rrni \\ ^ s 
a notation for the binomial coefficient which occupies less space in Latex "equation mode" than 
the classical notation 



(2.23) 



Now, we state directly the final theorem, without further inductive or intuitive information. 
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Theorem 2.24. For k > 1, we have: 



(2.25) 



K+l 



Y. = y.- + E E E E 

c£=l 1<Ai<---<A c ;<k- /ii>l,...,/i^>l /y.i AiH h/J-d^d^^+l 

«;•••(«- jtixAi i-id^d + 1) 



} , T"-Clil- - 



-^d^d N*d 



«•••(«- /UiAi ^<jA d + 2)0uiAi H h /x d Ad) 



(Ai!)^ Mi!---(Ad!)^ Md' 



« — /* l A i - 



^d A d + 1 yMlH l-Mrf- 1 



Once the correct theorem is formulated, its proof follows by accessible induction arguments 
which will not be developed here. It is better to continue through and to examine thorougly the 
case of several variables, since it will help us considerably to explain how we discovered the 
exact values of the integer coefficients ^ /il ' Xl) -"' (/i '' ,A * ) and 



2.26. Verification and application. Before proceeding further, let us rapidly verify that the 
above general formula (2.25) is correct by inspecting two instances extracted from Y5. 

Firstly, the coefficient of (yi) 3 ys in Y5 is obtained by putting k = 5, d = 2, Ai = 1, fi\ = 3, 
A 2 = 3 and /x 2 = 1 in the general formula (2.25), which yields: 



(2.27) 



n 5 ■ 4 ■ 3 ■ 2 • 1 ■ 6 

(l!) 3 3! (3!) 1 1! y3 



-20 X„. 



This value is the same as in the original formula (2.10): confirmation. 

Secondly, the coefficient of yi(y 2 ) 2 m Y 5 is obtained by k = 5, d = 2, Ai 
A 2 = 2 and /i 2 = 2 in the general formula (2.25), which yields: 



1, Mi = 1. 



(2.28) 



5-4-3-2-1 
(l!) 1 1! (2!) 2 2! 



5-4-3-2-5 
(l!) 1 1! (2!) 2 2! 



■<ij- 



75 X 



■I'.y- 



This value is the same as in the original formula (2.10); again: confirmation. 

Finally, applying our general formula (2.25), we deduce the value of Ye without having to 
use Y5 and the induction formulas (2.4), which shortens substantially the computations. For the 
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pleasure, we obtain: 

' y 6 = y x6 + [ey x5y - x xB ] yi + [i5^ v -6X x5y ] ( yi ) 2 + 



(2.29) 



20^3 



(yif + [isy X 2 yi 



■ 20 P^ x 3y3 



(y 



~6y xy5 - i5x x 2 yi ] (yi) 5 + [y yB -6X xy5 ] { yi f + [-x yS ] ( yi ) 7 + 

15 y x 4 y -6X x5 ] j/2 + [60y x s y 2 - 45 X x i y ] 2/12/2+ 

90y xV - 120^3^] ( yi fy2 + [<50y xy 4 - 150^3] (i/i) 3 t/ 2 + 

15^5 - 90^ 4 ] ( yi ) 4 y 2 + [-21^5] (yi) 5 y 2 + 

A5y x 2 y 2 - 60X x3y ] (y 2 ) 2 + [90y xy3 - 225X x2y2 ] yi(y 2 ) 2 + 

45 3V - 270 A^] (yi) 2 (y 2 ) 2 + [-210 A^*] (yi) 3 (y 2 ) 2 + 

lby y3 - 90X xy2 ] (y 2 f + [-105 ^3] yi( V 2) 3 + 

20y x3y - 15X xi ) y s + [60^2 y2 - 80^] 2/12/3+ 

60y xy3 - 150 A^] { yi ) 2 y 3 + [20 3V - 120X xy3 ] (yi) 3 y a + 

-35X yi ] (yi) 4 j/3 + [60y xy 2 - 150X x 2 y ) 2/22/3+ 

60 3V - 360 A"^] 2/12/22/3 + [-210 AV] (yi)' 2 y2y 3 + 

105 AV] (y 2 ) 2 y s + [wy y 2 - 60 X xy ] (y 3 ) 2 + 
-70X y2 ] 2/1 (y,f + [I5y x 2 y - 20 A>] y 4 + 
30y xy 2 - 75X x 2 y ) j/12/4 + [15 3V ~ ^ X xy2 ] (yi) 2 y4+ 
-35 AV] (yi) 3 y4+ [15 3V -90X xy ] 2/22/4+ 
-105 X y2 ] yi2/22/4 + [-35 Ai] 2/32/4 + [63V, - 15 AV] Vs+ 
6 y y 2 - 36 X xy ] yi y 5 + [-21 X y2 ] (yi) 2 y 5 + [-21 X y ] 2/22/5 + 
y y - 6 X y ] 2/c + [-7X V ] 2/12/6- 



2.30. Deduction of the classical Faa di Bruno formula. Let x, y G K and let 3 = g{x), 

f = f(y) be two C°°-smooth functions K — > K. Consider the composition h := / o g, namely 
h(x) = f(g(x)). For A 6 N with A > 1, simply denote by g\ the A-th derivative and 

similarly for h\. Also, abbreviate f\ := ^pr- 

By the classical formula for the derivative of a composite function, we have h\ = f\g\. 
Further computations provide the following list of subsequent derivatives of h: 



(2.31) 



hi = 
h 2 = 
h 3 = 
hi = 
< h 5 = 



/2 (.91 ) 2 + A .92, 
h (5i) 3 + 3/ 2 5i32 + /i33, 

A (si) 4 + 6/3 (si) 2 32 + 3/2 (g 2 ) 2 + 4 / 2 ffi 33 + /1 34, 
/s (3i) 5 + 10 / 4 (3i) 3 32 + 15 h (3i) 2 .93 + 10 / 3 31 (.9 2 ) 2 + 
+ 10 h 32 .93 + 5 f 2 31 .94 + /1 35, 

f 6 (31 ) 6 + 15/5 (3i) 4 32 + 45/ 4 ( 5l ) 2 (3 2 ) 2 + 15/3 (32) 3 + 
+ 20 U (3i) 3 33 + 60 h 3i .92 33 + 10 f 2 C93) 2 + 15 fa {gif 34+ 
+ 15 f 2 32 34 + 6 f 2 31 35 + A 3 6 . 

Theorem 2.32. For every integer k > 1, f/ze K-f/i derivative of the composite function h = f o 3 
may fee expressed as an explicit polynomial in the partial derivatives of f and of g having integer 
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coefficients: 



(2.33) 



d K h 
dx K 


K 

= E 


E 


E E 






l<Ai <-..<Ad<K 


AH>l,---^ici>l AtiAiH h/i t jA 1 j=« 








k! 








(Ax 









This is the classical Faa di Bruno formula. Interestingly, we observe that this formula is 
included as subpart of the general formula for Y K , after a suitable translation. Indeed, in the 
formulas for Yi, Y2, Y3, Y4, Y5, Y6 and in the general sum for Y K , pick only the terms for 
which |Ui Ai + • • • + jUdAd = K and drop X, which yields: 



(2.34) 



EE E E 

d=l l<X 1 <-<X d <K fix>l,...,fj,d>l \-^ d X d =K 



The similarity between the two formulas (2.33) and (2.34) is now clearly visible. 

The Faa di Bruno formula may be established by means of substitutions of power series 
( IF 19691 . p. 222), by means of the umbral calculus ( ICS 19 96 1), or by means of some induc- 
tion formulas, which we write for completeness. Define the differential operators 

d ( d \ 

d d ( . d . d \ 

Fs ■= 92 5 h gz 5 1-51/2 + h -^r , 

(2.35) dgi dg 2 \ dfi df 2 J 



d d 

Fx := gi V g- A h ■ 

dgi agi 



d , /. d d 
9x -5 K 31 [h -wj- + h 



+ h 



dfx 



Then we have 



(2.36) 



hi = F 2 (hi), 
h 3 =F 3 (h 2 ), 



h x =F x (h x _ 1 ). 



§3. Several independent variables and one dependent variable 

3.1. Simplified adapted notations. As announced after the statement of Theorem 2.24, it is only 
after we have treated the case of several independent variables that we will understand perfectly 
the general formula (2.25), valid in the case of one independent variable and one dependent 
variable. We will discover massive formal computations, exciting our computational intuition. 

Thus, assume n > 1 and m = 1, let k £ N with k > 1 and simply denote (instead of (1.2)) 
the jet variables by: 

(3-2) (x\y,yi 11 yi 1 _i 2 , . . . iVi-i^,. ■ 
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Also, instead of (1.30), denote the «-th prolongation of a vector field by: 

n 

- E Y 4l , 

i 

d 



(3.3) i 



il r\ £5 £1 

' dx l dy ^—^ 11 dvi, 

i—l ii — 1 ii,Z2 — 1 



' 2 %l,i 2 



E Y - 



»i,ta, ...,i„=l 



The induction formulas are 
Y 4 , : 



(3.4) 



= D i(y)-J2 D l( xk 



Vk, 



k=l 



Y ilti2 :=D 2 i2 (Y il )-^Dl 2 {X k 



Vil ,fc; 



k=l 



IK (Y .... ) X D l Wi 



fe=l 



we have 



where the total differentiation operators are defined as in (1.22), dropping the sums 2™=i 
and the indices j'. 

3.5. Two instructing explicit computations. To begin with, let us compute Y» I . With = 
have: 

n 

y u = Ai (y) - E ^ C** 1 ) 

fci=i 

n n 

fei=i fci=i 

Searching for formal harmony and for coherence with the formula (2.6) i, we must include the 
term y y inside the sum J2k 1= i H Using the Kronecker symbol, we may write: 



(3.6) 



(3-7) y y y H = [%yv\V*i- 

fei=i 

Also, we may rewrite the last term of (3.6) with a double sum: 

n n 

(3-8) - Xy 1 Vii Wfe = E WiMfe" 

fcl = l fcx,fe2 = l 

From now on and up to equation (3.39), we shall abbreviate any sum Ylk=i f rom 1 to n as ^2 k . 
Putting everything together, we get the final desired perfect expression of Yj x : 



(3.9) 



Yn = + E K yy x %] + E K 1 ^ 



fcl ,&2 



This completes the first explicit computation. 

The second one is about Y^ ^ 2 . It becomes more delicate, because several algebraic transfor- 
mations must be achieved until the final satisfying formula is obtained. Our goal is to present 
each step very carefully, explaining every tiny detail. Without such a care, it would be im- 
possible to claim that some of our subsequent computations, for which we will not provide the 
intermediate steps, may be redone and verified. Consequently, we will expose our rules of formal 
computation thoroughly. 
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Replacing the value of Yi just obtained in the induction formula (3.4)2 and developing, we 
may conduct the very first steps of the computation: 



/ UiiM 

fci 

3a* +Vi2 dy + ^ V ' 2M dy k 

fci 



(3.10) 

( d 



V dx^ 



-j (y*n +E [€ i yy-^X]y^+ 
+ E [- 5 n x v 2 ] ) - E [ x % + ^ < x ] 

fcl,fc2 / fcl 

I + E K 1 yy - X %] V*r + E [~ S Z X v 2 } V-rV» I + 

\ fcl kx,k-2 J 

\y X H + E K 1 yy - x %] + E < 2 ] ) + 

\ fci fc'i,fc2 ) 



d_ 

dy 



+ (? ^ ^ 



- ] f^n+E % ~ ] + E X v 2 } W>rM» ) + 

'7 \ fci fci,fc 2 / 

+ E [~ X %] ykl ^ + E f - ^ 1 ] VnVnM 

fcl &1,&2 

+ ^'ijfe + E [^i 1 y w* ~ f^i^a + E [~ s il x yy) Vk±yk 2 yi2+ 

fcl ^1^2 

+ E K 1 yy ~ X X] V^M + E X v 2 ] V**V**M + E X v 2 ] V»iV*M+ 

fcl fc'l ,fc2 

+ EK4K* 1 +E[-< 1 ] 



fcl fc'i ,fc2 fcl jfc2 



fcl fc i 

Some explanations are needed about the computation of the last two terms of line 11, i.e. about 
the passage from line 7 of (3.10) just above to line 11. We have to compute: 



(3.11) 

This term is of the form 



fcl Vkl / \ki,k 2 



(3-12) E A kiirM(E [ B fci,fc2]ffci2/fc 2 

\fci Wm 

where the terms B kl k 2 aie independent of the pure first jet variables y x k . By the rule of Leibniz 
for the differentiation of a product, we may write 

|E Afc i^| ( E [ B ki,k 2 ] VkiVk 2 



fcl / \ fel »&2 



(3.13) 



E [ B fci,fc 2 ]s/fc 2 E A k' ± a^-(ffci) + E [^.Jkalv*! (E A 
fei.fea \ fc; y M / fci.fea V fc 2 



' q (Uk 2 ) 



E [ B kl,k 2 ] Vk 2 A kl + i B kl,k 2 ] Vkt A k2 - 

k\,k 2 k 1 ,k 2 
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This is how we have written line 1 1 of (3.10). 

Next, the first term y x i ly yi 2 in line 10 of (3.10) is not in a suitable shape. For reasons of 
harmony and coherence, we must insert it inside a sum of the form Ylki H ^ fc i ■ Hence, using the 
Kronecker symbol, we transform: 



(3.14) 



fel 



Also, we must "summify" the seven other terms, remaining in lines 10, 11 and 12 of (3.10). 
Sometimes, we use the symmetry y^ki = 2/fci,i 2 without mention. Similarly, we get: 



fei 

E 

fel ,&2 

£[<^ 



fei .fe2 



A': 



yk-LVk^y i2 



e 



fel ,^3 



V K 1 

W ~ 'Vi 



2/fei, 



yk ± yk 2 
yk 2 yk 3 , 



fei ,fe 2 



(3.15) 



E R 1 ^ 



yk 2 yk\,i 2 



k\ ,k 2 



E hM 1 



fel ,fe2 



E -<t« 



y^i yk 2 ,i2 

^3 V^l 



yk 1 yk 2 ,k 3 i 



E « 
£[-<*£ 



yk\yk 2 ,i 2 



yki,h 



E 



<M 2 



Aii ,Ai2 ,fes 



yk t yk 2 Mi 



E 



fel >fe2 



yfei ,fe2 j 



E t"^ 1 ] = E t^^y 2 ] ^J/fe 2l ii 



fei 



E -<t« 



2/fei2/fe 2 ,fe 3 - 



fel,fe2,fe3 



In the sequel, for products of Kronecker symbols, it will be convenient to adopt the following 
self-evident contracted notation: 



(3.16) 



S? 1 Sf 2 = 5? 1 '? 2 ; generally : Sf 1 5f 2 ■ ■ ■ Sf x = 5, 1 



12, ■■■ ,l\ ' 



16 



JOEL MERKER 



Re-inserting plainly these eight summified terms (3.14), (3.15) in the last expression (3.10) of 
Yi lt i 2 (lines 10, 11 and 12), we get: 



ki fcl>*:2 



+ E K 1 ^ J vm + E [<* ,t y w - € 2 X % y ] Wi + 
fc i m fc i> fe 2 m 







(3.17) 



+ E HiV**. 



-0 



-0 



ki,k 2 ,k 3 



yy VkiVk 2 Vk 3 



+ 



fei ,k2 



fcl,fc 2 ,fc 3 



+ 



E [-6%x£]v** + 





ki ,k 2 



^0 

E [-<y£ 3 ** 2 K 

ki ,^2,^3 



-0 



-0 



Next, we gather the underlined terms, ordering them according to their number. This yields 6 
collections of sums of monomials in the pure jet variables: 

Yu.ia = + E [ S il y ^y + 5 iz y^y ~ X % x i 2 \ W>i + 

fei 

+ E [ 5 tt yyy - # *X - **% X X] yk ^+ 

fei 

(3.18) + E [-^iVfr ^m] ykiVk 2 Vk 3 + 

k 1 ,k 2 ,k 3 

+ E KiV&^- a &^-< 2 ^4]^+ 

fci,fe2 
ki,k 2 ,k 3 

To attain the real perfect harmony, this last expression has still to be worked out a little bit. 
Lemma 3.19. The final expression of Yi t i 2 is as follows: 



(3.20) 



Yi 



fci 

cfci x>fc 2 JC&1 



V k fel f 2 y w - (S* 1 — sf 1 x% 



Vk x yk 2 + 



f fci,fc 2 ^fc 3 



2/fci2/fc 2 yfc 3 - 



fcl ,&2 ,fc3 



fcl jfe2 



fcl,fc2 vfc3 ffc3,fcl 1^2 ffc2,fc3 -yfcl 



fcl ^2 jrCS 



Vk l Vk 2 ,k 3 - 



Proof. As promised, we explain every tiny detail. 

The first lines of (3.18) and of (3.20) are exactly the same. For the transformations of terms 
in the second, in the third and in the fourth lines, we use the following device. Let Tki,k 2 be 
an indexed quantity which is symmetric: Tk 1 ,k 2 = Yfc 2 ,fci- Let Ak u k 2 be an arbitrary indexed 
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quantity. Then obviously: 

(3-21) ^ -<4*i,fc2 Tfc ljfe2 = 4 2 ii Tfcj.tj. 

^1,^2 k\,k 2 

Similar relations hold with a quantity ^ ; A which is symmetric with respect to its A indices. 

Consequently, in the second, in the third and in the fourth lines of (3.18), we may permute freely 
certain indices in some of the terms inside the brackets. This yields the passage from lines 2, 3 
and 4 of (3. 18) to lines 2, 3 and 4 of (3.20). 

It remains to explain how we pass from the fifth (last) line of (3.18) to the fifth (last) line 
of (3.20). The bracket in the fifth line of (3.18) contains three terms: [-T x - T 2 - T 3 ]. The term 
T3 involves the product S^'^, which we rewrite as 6^'^, in order that i\ appears before £2. 
Then, we rewrite the three terms in the new order [— X2 — T3 — T{\, which yields: 



(3-22) [-%X'x$*-*%$x!r-%;%x* 1 

^1,^2,^3 



It remains to observe that we can permute &2 and k$ in the first term — T2, which yields the last 
line of (3.20). The detailed proof is complete. □ 

3.23. Final perfect expression of Y il i2 i3 . Thanks to similar (longer) computations, we have 
obtained an expression of Y il i2 i3 which we consider to be in final harmonious shape. Without 
copying the intermediate steps, let us write down the result. The comments which are necessary 
to read it and to interpret it start just below. 

Yi 1 ,i 2 , 13 = J'l'll'Jl'S + ^>ii y x i 2x i 3y + x l 3 y + 3^*1 a *2 v — ^j'li;'!^] 2/^1 + 

fcl 

_l_ \ " [rfcl>fe2 -y . I Sjk 1 ,k 2 -y , Sjki,k 2 -y _ 

' / j ii,*2 x l 3 y% I 13 , i\ *s x l 2 y 2 ' i 2 1 *3 **x l ly 2 



(3.24) 



ki ,k 2 

~' J il "~x'2x i 3y "*2 'Vll'S» "*S 'Vli'Sj 
4- [AA-'ji; AAA _ifci,fc2vfc3 _sfci,fc2v fe 3 1 1,1,4- 

^1 > fe 2i fc 3 

+ ["^iV'^'is 3 -^v] ykiVk 2 Vk 3 yk i + 

ki,k2,k3,k^ 

+ V* |"(5 fel ' fe2 ' fe3 V 2 + 5 fc3 ' fel ' fe2 v 2 + s k2 ' k3 ' kl y 2 - 

fcl,/C2,/c3 



x'3y 


_ £^2^3 


x'2y 


- <£ 2 i 3 

'1 j b 3 x y 


x'ly 


_ £k 2 ,k 3 p^kx 

*2> *3 x i ly 



J Vk 1 yk 2 M + 



1 \ ^^1,^2^3 V fc 4 £k 2 ,k 3 ,ki yk.4 £k^,k 2 ,k 1 yk 4 

+ 2^ [ °ii,i2,<a y 2 ~°iii<a,<a y 2 _0 ii,*2^3 *V ~ 

fc l: fc 2 ? fc 3' fc 4 

_ J /= 3 ,fc 4 ,fc 1 x ki _ gkaMM x h 2 _ d kuk 3 ,k 4 x k^ 



ki,k 2 ,ks,k4 
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kl,k 2 ,k 3 

, \ ' S k l, k 2 , k 3 V k 4 Si k i, k l, k 1 V fc 3 : k 3,k i ,k 1 -yk 2 X k 2< k 3' k 4 V k l \ „, „, 

k l, k 2, k 3 , k 4 

3.25. Comments, analysis and induction. First of all, by comparing this expression of Y^ ^ ^ 
with the expression (2.8) of Y3, we easily guess a part of the (inductional) dictionary beween 
the cases n = 1 and the case 77. > 1. For instance, the three monomials [-](yi) 3 , [•] 7/12/2 and 
[•] (yi) 2 7/2 in Y3 are replaced in Y^ ^.i, by the following three sums: 

(3.26) ]T [^ykiVksVka, Yl Hffci2/fc2,fc 3 > and E [■]ykiVk 2 yk 3 ,k i . 

Similar formal correspondences may be observed for all the monomials of Yi, Y^, of Y2, 
Yi 1; i 2 and of Y3, Yi 1) i 2) i 3 . Generally and inductively speaking, the monomial 

0.27) {■}(yx 1 T 1 ---(yx d f d 

appearing in the expression (2.25) of Y K should be replaced by a certain multiple sum generaliz- 
ing (3.26). However, it is necessary to think, to pause and to search for an appropriate formalism 
before writing down the desired multiple sum. 

The jet variable y\ 1 should be replaced by a jet variable corresponding to a Ai-th partial deriv- 
ative, say j/fcj fc A , where fc 1; . . . , k\ 1 = 1, . . . , n. For the moment, to simplify the discussion, 
we leave out the presence of a sum of the form kx . The /ii-th power (y\ 1 ) Ml should be 

replaced not by (yki,...,k\ , but by aproductof fxi different jet variables y kl k Xl oflength 

A, with all indices k a = 1, . . . , n being distinct. This rule may be confirmed by inspecting the 
expressions of Y ii; of Y ilji2 and of Yj 1 ,, 2 ,j 3 . So yk 1 ,...,k Xl should be developed as a product of 
the form 

(3-28) yki,...,k\ 1 Vk\ 1 +i,—,k2\ 1 yk^ t i 1 -i)x 1 +i,--;k t , 1 \ 1 > 

where 

(3.29) ki,...,kx 1 ,...,k fil Xi =l,...,n. 

Consider now the product (yx^ 1 {yx 2 Y 2 ■ How should it develope in the case of several inde- 
pendent variables? For instance, in the expression of » 2 » 3 , we have developed the product 
(2/1) 2 J/2 as yk 1 yk 2 yk 3 ,k i - Thus, a reasonable proposal of formalism would be that the product 
(yxi)^ 1 (2/a 2 ) A ' 2 should be developed as a product of the form 

(3 30) ffcij— ,feAi ykx 1 +i,-;k 2 x 1 2/fe( M1 -i)Ai+i!— l^lAi 

yk l i 1 \ 1 + l,---,k l i 1 \ 1 - l -\ 2 yk tll \ 1 +( f j, 2 -i)x 2 +i,---,k lll \ 1+l i 2 \ 2 i 

where 

(3.31) fei, . . . , kx t , . . • , k lll x l , . . • , k^ 1 x 1 +fi 2 X2 = 1) ■ ■ ■ > n - 

However, when trying to write down the development of the general monomial 
(j/Ai) Ml {yx- 2 Y 2 ■ ■ ■ (yx d ) fJ ' d , we would obtain the complicated product 

Vki,...,kx 1 Vk\ 1+1 ,...,k 3 x I '" yk( fll -i)\ 1 + i,---, k ^ 1 \ 1 

V k M X 1 + l>---> k n X 1 +X 2 ■ ■ ■ 2/fen 1 A 1 + (n 2 -l)A 2 + l.---. fc niA 1 + t i2^2 

(3.32) 



y k i+ 1 \ 1 +--+v. d _ 1 \ d ^ 1 +i>---> k v 1 \ 1 +--+^ d _ 1 \ d _ 1 +\ d 

y k f i 1 \ 1 +-+ f i d _ 1 \ d _ 1 +( Md -i)\ d +i,---> k v 1 \ 1 +--+v- d \ d ■ 

Essentially, this product is still readable. However, in it, some of the integers k a have a too 
long index a, often involving a sum. Such a length of a would be very inconvenient in writing 

k k 

down and in reading the general Kronecker symbols S^ 1 which should appear in the final 
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expression of Yj lt . One should read in advance Theorem 3.73 below to observe the presence 
of such multiple Kronecker symbols. Consequently, for a = 1, . . . , /iiAi, . . . , /xiAi + • • • + 
lid^d, we have to denote the indices k a differently. 

Notational convention 3.33. We denote d collection of pd groups of Ad (a priori distinct) integers 
k a = l,...,nby 

^1:1:1) ■ ■ ■ j &l:l:Ai ; ■ ■ ■ 3 ^l:/-ti:l j ■ ■ ■ 3 ^l:/Xi:Ai ) 

" * ' > v ' 

Ai Ai 



/'1 

fe:l:l, ■ ■ ■ , fe2:l:A 2 ) • ■ • ) &2:^ 2 :1j ■ • • ) k2:fj, 2 :\ 2 



(3.34) 



A 2 



/<■-> 



kd:l:l, ■ ■ • , kd:l:\ d , ■ ■ ■ , &d:^ d :l, ■ ■ ■ , fcd:/j d :A d ■ 



/'./ 

Correspondingly, we identify the set 

(3.35) {1, . . . , Ai, . . . , /LtiAi, , jUiAi + ^ 2 A 2 , , /LtiAi + a* 2 A 2 H + Pd^d} 

of all integers a from 1 to p\X\ + p-2^2 + • • • + Pd^d with the following specific set 

(3.36) {1:1:1, . . . , l:l:Ai, . . . , hpv.Xi, . . . , 2 : p 2 : A 2 , . . . , d : fid ■ Ad}, 

Xl 



M1 A 1+M2>>2 



MlXi+M2->>2H hMejA,j 

written in a lexicographic way which emphasizes clearly the subdivision in d collections of pd 
groups ofXd integers. 

With this notation at hand, we see that the development, in several independent variables, of 
the general monomial (yx^ 1 ■ ■ ■ {yx d Y d may be written as follows: 

(3.37) **' 2/^1:^1 :1)"'J^1:M1 '"' V^d: 1:1 )•") ^d: 1: X d V^d : p d : 1 j ■ ■ ■ , &d: p d : A d ' 

Formally speaking, this expression is better than (3.32). Using product symbols, we may even 
write it under the slightly more compact form 

(3-38) \\ Vk l! u 1 ,l,...M: V1 :X 1 il Vk*» d :L-,**.~ d * d - 

l<Vl<fj,i l<i> d <n d 

Now that we have translated the monomial, we may add all the summation symbols: the 
general expression of Y K (which generalizes our three previous examples (3.26)) will be of the 
form: 
(3.39) 

y k = ^w. + e E E E 

d— 1 l<Ai<---<A d <K (Ui>l,...,)Lid>l Mi^iH r-MdA d <K,-f-l 

n n n n 

E •■■ E E - E 

m n 

l<I/l</il l<V d <[J, d 
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From now on, up to the end of the article, to be very precise, we will restitute the bounds J2k=i 
of all the previously abbreviated sums J2k- This is justified by the fact that, since we shall deal in 
Section 5 below simultaneously with several independent variables (x 1 , . . . , x n ) and with several 
dependent variables (y 1 , . . . , y m ), we shall encounter sums Y^iLx> not to be confused with sums 

En 
k=l- 

3.40. Combinatorics of the Kronecker symbols. Our next task is to determine what appears 
inside the brackets [?] of the above equation. We will treat this rather delicate question very 
progressively. Inductively, we have to guess how we may pass from the bracketed term of (2.25), 
namely from 



k - ■ ■ (k — /ilAi 



Hd\d + 1) 



~l i d x d i/ClH - 



(3.41) 



(AiO^/i!!---^!)^! x ddy 

K ■ ■ ■ (k - fiiXi Hd^d + 2)(/iiAi H h^dAd) 



(Ai!)^ pi!--. (Ad!)** n d \ 



• X~ 



-f d x d + 1 d 



to the corresponding (still unknown) bracketed term [?]. 

First of all, we examine the following term, extracted from the complete expression of 



Y iui2 , i3 (first line of (3.24)): 

71 



E 

fei=i 



x z 2x t 3y 



S kl y h 



x l lx l 3 y 



5 ki y i, 



x z 2y 



x l i x l 2 x l 3 



2/fei- 



Here, the coefficient [3 34 2 j/ ~ ^x 3 ] of the monomial j/i in Y3 is replaced by the above brack- 
eted terms. 

Let us precisely analyze the combinatorics. Here, X x a is replaced by X i2 i3 , where the 
lower indices ii, i2, 13 come from Y^ ^ ^ and where the upper index fci is the summation index. 
Also, the integer 3 in 3 34 2 j, is replaced by a sum of exactly three terms, each involving a single 
Kronecker symbol 6\, in which the lower index is always an index i = i\, £2,23 and in which 
the upper index is always equal to the summation index k\. By the way, more generally, we 
immediately observe that all the successive positive integers 



(3.43) 



1,3,1,3,3,1,3,1,3,3,3,9,6,3,1,3,4 



appearing in the formula (2.8) for Y3 are replaced, in the formula (3.24) for Y^ ^ ^, by sums 
of exactly the same number of terms involving Kronecker symbols. This observation will be a 
precious guide. Finally, in the symbol S kl , if i is chosen among the set {i\,i2, 13}, for instance 
if i = i\, it follows that the development of 34 2 y necessarily involves the remaining indices, for 
instance 34^ ^3 y . Since there are three choices for i = ii, t2, «3, we recover the number 3. 
Next, comparing [y yy — 2 X xy ] (j/1) 2 with the term 



(3.44) 



n 

E [* 



sf 1 x\ 



sf 1 x\ 

22 x l iy 



fel ,fc 2 



extracted from the complete expression of i 2 (second line of (3.18)), we learn and we guess 
that the number of Kronecker symbols before y x ~t y & must be equal to the number of indices k a 
minus 7. This rule is confirmed by examining the term (second and third line of (3.24)) 



(3.45) 



st^y^yi + s^yx' 



X l 2 y 2 



5 klM y <, 



k\ ,fc 2 



_ rfci x k 2 



x l 2 x l 3 y 

developing [3y xy 2 -3X x 2 y ] ( yi ) 2 . 



S k! X k 



X 1 ! 3^3 y 



§ ki X k 



x z i x l 2 y 
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Also, we may examine the following term 

n 



ki,k 2 — 1 



(3.46) 



, rki,k 2 i r«l,K2 . _l_A Kl ' K2 "W . 

~ l ~ *2,l3 "r°i 2 ,i 4 ->Vlx*3y2 "r i3,i 4 •'s'li 



cki.k 



kl,k 2 



3*1 3; l 2a; l 3^ l 4y Z2 a; ! l a;^ j; ! 3 y ^3 a^ljE^xMy 



*4 a; 1 X 2 X 3 y 



extracted from Y^ ^ ^ ^ and developing [6y x 2 y 2 — ^^x s y] (yi) 2 - We would like to mention 
that we have not written the complete expression of Y^^^^, because it would cover two and 
a half printed pages. 

By inspecting the way how the indices are permuted in the multiple Kronecker symbols of 
the first two lines of this expression (3.46), we observe that the six terms correspond exactly to 
the six possible choices of two complementary ordered couples of integers in the set {1, 2, 3, 4}, 
namely 

{1,2}U{3,4}, {1,3}U{2,4}, {1,4}U{2,3} > 

{2, 3} U {1,4}, {2,4} U {1,3}, {3, 4} U {1,2}. 

At this point, we start to devise the general combinatorics. Before proceeding further, we need 
some notation. 

3.48. Permutation groups. For every p E N with p > 1, we denote by & p the full permutation 
group of the set {1, 2, . . . ,p — l,p}. Its cardinal equals pi. The letters a and r will be used to 
denote an element of & p . If p > 2, and if q € N satisfies 1 < q < p — 1, we denote by 6^ the 
subset of permutations a e 6 p satisfying the two collections of inequalities 

(3.49) a(l) < <i{2) < ■ ■ ■ < a{q) and a{q + 1) < a{q + 2) < • • • < a{p). 

The cardinal of S« equals C« = g , £_ q y . 

Lemma 3.50. For k > 1, the development of (2.20) to several independent variables 
(x 1 , . . . , x n ) is: 



~Yi lt i 2 ,...,i K — y x i-L X H--x*K + El 

71 

+ E 

fc l >&2 — 
71 

+ E 



a; ■ ■ -x ' 



^1 !^3 — 1 



(3.51) 



661 

V"* £ fc l> k 2 -y . _ rfcl V fe 2 

V(1)>V(2) J V T C3) ...a,*T(K) v a Ml) x i r(2)... x i T 

5| rSS; 

V"* ifcl, k 2 , k 3 y 
V(l),V(a).'r(!) 

66 J 

\ ^ I^li ^2 ^^3 
Z^ Ml),V (2 ) a! iT (3)... a! i TW 1/ 



3/fci + 



Vk 1 Vk 2 + 



t£6 



Vk 1 yk 2 yk 3 + 



■ + 



+ E 

fe 1 ,...,*: re =i 



On,..., !k JV - Z. d i Tt i).-.M»-i) 



res; 



Vki ■ ■■Vk K + 



n 

+ E [~ S il'Z'in X y£ +1 ] VH'-' Vk K Vk K+1 + remainder. 



fci , . . . ,fc K , k K —1 
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Here, the term remainder collects all remaining monomials in the pure jet variables yu x , 



3.52. Continuation. Thus, we have devised how the part of ^ which involves only the jet 

variables yk a must be written. To proceed further, we shall examine the following term, extracted 



from Y,, 



(3.53) 



(lines 12 and 13 of (3.24)) 



E 

fci ,k->,k<>,ki 



cki,k2,ks -yki 
11,12,13 y 2 



_ck3,k 4 ,ki yk 2 
-°ii,i 2 ,i 3 y 2 



ck2,k3,ki -yki 

41,42,13 y 2 



jfe3,fel,fc4 -^k 2 



ck3,k2,ki -yki 

81,42,43 y 2 



j/ci,fc3,fc4 ^>&2 



VkiVk2Vk3M, 



which developes the term [—6 X y i\ (yi) 2 y2 of Y3 (third line of (2.8)). During the computation 
which led us to the final expression (3.24), we organized the formula in order that, in the six 
Kronecker symbols, the lower indices 11,1-2,13 are all written in the same order. But then, what 
is the rule for the appearance of the four upper indices ki, fe, kz, ktP. 

In April 2001, we discovered the rule by inspecting both (3.53) and the following complicated 
term, extracted from the complete expression of Yi Il i 2i i 3i j 4 written in one of our manuscripts: 

E " 



(3.54) 



^ki,k2,k3 -vi 
il, 42, 13 ^x'iy 2 


+ 


rk 2 ,ki,k3 -y 
ii, 42, 43 Jx % iy- 


+ <5 


k2,k3,ki -vi , 

41, 42, 43 ^l' 4 !/ 2 " t 




, ckl,k 2 
' ii, 4 2 


k . 3 y , 

44 •Sx 


3y 2 


1 $k2,ki,k3 -vi 

' 41 , l 2 , 44 <SX 


*3y2 


, fikzMM y , 

^ 4 2 , 44 • , r'r T 


, rki,k 2 
' ii, 43 


k3 y , 

44 •sx' 


2y2 


1 T^'2 ,^'l ,^3 "\1 
' °4i , 43, 44 


*2y2 


1 xk 2 ,k 3 ,ki . 

+■ d ti, i 3 , 44 >X , 2 3 , 2 + 


, Sjki,k 2 

^ 42, 43 


H y , 

44 -fx* 


ly2 


1 T^'2 ,^'1 ,^3 "\1 
4 2 . i 3 , 14 


i ly 2 


1 ^^2,^3,^1 -\; 

°42, 43, 4 4 ^ n 


y 2 - 


cki,k 2 

il, 4 2 


x\ 


14 y 


rfc 2 ,fel vfc3 
41, 4 2 x i 3x 


iy 


srk 2 ,k3 -yki 

il, 42 l'31'lj 




ijki,k 2 

ii, 43 


X z 2 x 


14 y 


_ rfc 2 ,fci yk 3 
il, 43 i'2i 


4y 


Sjk2,k3 yki 
il, 43 l'2l'4j 




_ Sjki,k 2 

ii, 44 


A* ^ 3 

X z 2 x 


>3y 


— A fe2,fel Af fe3 

U 4i, 44 a;i2a; 


3 V ~ 


r&2,&3 'yki 

il, ii x i 2x i 3y 




<rki,k 2 
%, 4 3 


x\ 

X 1 X 


14 y 


42, 43 X l lX 


4y 


Sf 2 ' k3 X k z \ u 

42, 43 x z ix z iy 




_ <;ki,k 2 

12, 4 4 


x k i 

X 1 X 


'3y 


rfc 2 ,fci -v>fc 3 
4 2 , 4 4 ^1*1! 




Sjk 2 ,k 3 -yki 

42, 44 x i lx i 3y 




_ <rk 1 ,k 2 
43, i 4 


x k i , 

ari x l 


2 11 


_ ^k 2 ,k 1 x>k3 

43, 4 4 iHf 


2 y ~ 


£^2,^3 V^l 
43, 44 x i lX i 2y 


Vki 



This sum developes the term [12 y xy 2 — 18 X x 2 y \ y\y^ of Y3 (third line of (2.9)). Let us explain 
what are the formal rules. 

In the bracketed terms of (3.53), there are no permutation of the indices £1, £3, but there is a 
certain unknown subset of all the permutations of the four indices ki , k 2 , k 3 , /c 4 . In the bracketed 
terms of (3.54), two combinatorics are present: 

• there are some permutations of the indices i\, £2, £3, £4 and 

• there are some permutations of the indices fci , fc 2 , £3. 

Here, the permutations of the indices i\ , £2, £3, £4 are easily guessed, since they are the same as 
the permutations which were introduced in §3.48 above. Indeed, in the first four lines of (3.54), 
we see the four decompositions 

(3.55) {£1, £ 2 , £3} U {£4}, {£i,«2,£4}U{£ 3 }, {£i,£3,«4}U{£ 2 }, {£2, £3, U} U {£1}, 

of the set {£1, £2, £3, £4}, and in the last six lines of (3.54), we see the six decompositions 

{£i,£ 2 } U {£3, £4}, U {£2, £4}, {£1,44} U {£2, £3}, 

{£2, £3} U {£1, £ 4 }, {£2, £4} U {£i,£3}, {£3, £4} U {£i,£ 2 }, 
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so that (3.54) may be written under the form 
(3.57) 



E 



V V X fc T(l)' fc x(2). fc x(3) y . _ : k T(l)' k T{2) V k T{3) 

Z^ Z^ V(1)>V(2)>V(3) y 2 Zj Z^ V(1)'V(2) x i T<,3) x i T{A) 



TgS? 



yk x yk 2 ,k 3 



where in the two above sums J2a-el' tne l etter a denotes a permutation of the set {1, 2, 3} and 
where the sign ? refers to two (still unknown) subset of the full permutation group 63. The only 
remaining question is to determine how the indices k a are permuted in (3.53) and in (3.54). 

The answer may be guessed by looking at the permutations of the set {k\, £4} which 

stabilize the monomial yk 1 yk 2 yk 3 ,k 4 in (3.53): we clearly have the following four symmetry 
relations between monomials: 



(3.58) 



yk 1 yk 2 yk 3 ,k 4 = yk 2 yk 1 yk 3 ,k i = yk x yk 2 yk 4 ,k 3 = y^ytuVkiMi 



and nothing more. Then the number 6 of bracketed terms in (3.53) is exactly equal to the cardinal 
24 = 4! of the full permutation group of the set {fci, k%, k^, k^} divided by the number 4 of these 
symmetry relations. The set of permutations a of {1, 2, 3, 4} satisfying these symmetry relations 



(3.59) 



2/fco-(l) 2/fecr(2) 2^(3) ,fccr(4) ~~ ^ fc 1 ^ fc 2 2/&3 , &4 



(2 1) (1 2) 

consitutes a subgroup of 64 which we will denote by #4 ' ' . Furthermore, the coset 



(3.60) 

possesses the six representatives 



(24), (1,2) . = g 4 ^(2,l),(l,2) 



(3.61) 



12 3 4 

12 3 4 

12 3 4 

3 4 12 



12 3 4 

2 3 14 

12 3 4 

3 14 2 



which exactly appear as the permutations of the upper indices of our example (3.53). Of course, 

(2 1^ fi 2^ 

the question arises whether the choice of such six representatives in the quotient ©4/^4 ' 
is legitimate. 

(2 1) (1 2) 

Fortunately, we observe that after conjugation by any permutation a € $)\' ' , we do not 
perturb any of the six terms of (3.53), for instance the third term of (3.53) is not perturbed, as 
shown by the following computation 



(3.62) 



[ r fc "(3). fc CT (2)>fc<r(l) yKw] „, „, „, _ 

[~ b il, «2, i 3 X V 2 \ VkiVk 2 yk 3 ,ki - 

k 2 ,k 3 , fc 4 

?k 3 ,k 2 ,ki y k(7 ( 4 )~\ 

~ Z^ l°n,'h, » 3 y 2 J ^-i(i) J/fc , T -i(2) J/ 'v-i(3)''v-i(4) 

kl,k 2 ,k 3 ,k 4 

sk 3 ,k 2 ,k-\ -t> fe cr(4)l 

= Z. [~ S h 'i2,il X y ' J VkxVhzVk 3 M 

k\ ,k 2 ,k 3 ,k 4 

thanks to the symmetry (3.59). Thus, as expected, the choice of 6 arbitrary representatives a G 
gf(2> 1 ).(i»2) j n thebracketec 
be written under the form: 



(2 1) (1 2) 

$ 4 ' ' in the bracketed terms of (3.53) is free. In conclusion, we have shown that (3.53) may 



(3.63) 



E 

k\ ,k 2 ,k 3 , k 4 



Er fe ff(l)> /i: t7(2) I Kcr(3) -^«<t(4) 



CT65 (2,D,(1,2) 



yk x yk 2 yk 3 M 



This rule is confirmed by inspecting (3.54) (as well as all the other terms of » 2 » 3 and of 
Yi lj i 2j i 3j i 4 ). Indeed, the permutations a of the set {fci, &2, ^3} which stabilize the monomial 
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Uk 1 yk 2 ,k 3 consist just of the identity permutation and the transposition of k% and ks- The coset 
©3/JQ3 1 ' has the three representatives 



(3.64) 



12 3 
12 3 



12 3 
2 1 3 



1 2 3 

2 3 1 



which appear in the upper index position of each of the ten lines of (3.54). It follows that (3.54) 
may be written under the form 



(3.65) 



E\ " rfe CT (l),fc CT (2).fe CT (3) y 

EST A~ fe < T < 1 >' fc < T < 2 > yk a (3) 

Z_^l V(1),V(2) x i T{3) x i T[i) 

aeei Te5 (i.D,(i,2) 



Vk l yk 2 M- 



3.66. General complete expression of Y< It . As in the incomplete expression (3.39) of 
Yi 1; .. , t i K , consider integers 1 < Ai < • • • < Xd < K and fix > 1,...,/^ > 1 satisfying 

fj,i\i + • • • + [idXd < « + 1. By ij^AH a > we denote the subgroup of permutations 

t € S^jXh M5„ * tnat l eave unchanged the general monomial (3.38), namely that satisfy 

,„ l<!/l</il l<V d <fl d 

(3.0/) 

1<I^1<^1 l<U d <fJ, d 

The structure of this group may be described as follows. For every e = 1, . . . , d, an arbitrary 
permutation er of the set 



(3.68) 



{e: 1:1, . . . , e: l:A e ,e: 2:1, . . . , e: 2:A e , ■ • • , e: /i e :l, . . . , e: ^ e :A e } 



which leaves unchanged the monomial 
(3-69) [] VK^,-,K 



r(e:^ e :l).---.' s 'iT(e:^e:^e 



1 J yk eMe:1 ,...,k e .. Ue 

l<v e <He 



l<v e <He 

uniquely decomposes as the composition of 

• [i e arbitrary permutations of the /i e groups of A e integers {e : v e : 1, . . . , e : v e : A e }, of 
total cardinal (A e !) Me ; 

• an arbitrary permutation between these \i e groups, of total cardinal /i e l. 

Consequently 



(3.70) Card 
Finally, define the coset 



(3.71) 
with 

(3.72) 



~(/ii,Ai),...,(/i<j,Ad) p. /f,\J*i,Ai),--->W;' 



(/ii,Ai),...,(/id,Ad) 



Card /^i,A 1 ),...,( M d,Ad)^ _ Card(6 AllAl+ ... +M dAd) 

I "/iiAiH YHd^d 



Card 55 



. (/Ji,Ai),...,(/jd,A d ) 
Mi^iH h^dAd 

(/X1A1 -I 1- /XdAd)! 

MiKAi!)" 1 •■•/i d !(A d !)^ 



EXPLICIT PROLONGATIONS OF INFINITESIMAL LIE SYMMETRIES 



25 



In conclusion, by means of this formalism, we may write down the complete generalization of 
Theorem 2.24 to several independent variables. 



Theorem 3.73. For every k > 1 and for every choice of k indices i\, 



in the set 



{1,2, 
(3.74) 



i}, the general expression of j K is as follows: 



Yii,...,^ —y x ii... x *„ + E y ' y_ E 

d=l l<Ai<..-<A d < K ni>l,...,M<i>l flAi + --+M<jAd<«+l 



E •■ E 

fc 1:1:1 ,...,fe 1:1:Al =l fcl: ( i 1 :l>'--,fcl: (il :A 1 =l 



E 



E 



~(MlAl) ( l i d ,\ d ) rg6 MlAl + -- + Hd A d 

"mi^iH \-v-d x d K 

gK(l:l:l)> — ' k a(l:lii:Xl)> — >Kld:li d :^ d ) 



E •■■ E 

k d:l:l,— <l"d:l:X d -l k d:n d :l k d^ d ->-d" 



qk-piXi Md^d+MiH V^dy 



E E 



, (Ml.*i),...,(n d ,A d ) T(? gA'lAl+-+|Jd A £i- 1 
^l^l + '-' + ^d^d 



7-(l:l:l)>-".K<r(l: M1 :A 1 ).---.«<r(<l:M c (:A d -l) 



QK- Md-^d + MlH YHd X k a{d:v. d ' > -d') 



, t( 1 )'"-' , t(*>1>0' — ''tCmi*I + — +(«d> < l-l) 9 2 ; i T(n 1 A 1 +-+ A1(i A ti ) . . . Qgf-T^) (py) M1 + " ' + Md ~ 1 
il "' il D k d:„ d :l,- ■ ; k d:v d :\ d - 



l<^d<Md 



3.75. Deduction of a multivariate Faa di Bruno formula. Let n e N with n > 1, let x = 
(x 1 , . . . , x") 6 K n , let g = g(x\ . . . , x") be a C°°-smooth function from K n to K, let y E K 
and let / = f(y) be a C°° function from K to K. The goal is to obtain an explicit formula for the 
partial derivatives of the composition h := f o g, namely h(x 1 , . . . , x") := f(g(x 1 1 . . . , x")). 
For A 6 N with A > 1 and for arbitrary indices ii, . . . ,i\ = 1, . . . , n, we shall abbreviate the 
partial derivative dx if. g dx i x by gi u ...,i x and similarly for hi u ...,i x . The derivative will be 
abbreviated by f\. 

Appying the chain rule, we may compute: 

= h [ff»i] , 

hii,i 2 ,i 3 = /3 [ffti 3*2 3*s] + /a [ffn 3*2, *3 + 3*2 3n,*3 + 9*3 3«i,i 2 ] + h [ffn.ta.ts] > 

/Hl,<a,<3,i4 ~ /4 [ffil 3*2 3*3 3*4] + /3 [3*2 3*3 3*1, *4 + 3*3 3*1 3*2, *4 + 3*1 3*2 3*3, *4 + 

+3*1 3*4 3*2, *3 3*2 3*4 3*1, *3 + 3*3 3*4 3*1, 12] + 
+ /2 [3*1 ,*2 3*3, *4 ~l" 3*1 ,*3 3*2, *4 "I" 3*1,14 3*2, '3] "I" 

+ /2 [3*1 3*2,*3,*4 "I" 3*2 3*1,*3,*4 "I" 3*3 3*1,*2,*4 "I" 3*4 3*1,*2,*3] + 

+ A [3*i 



(3.76) 



Introducing the derivations 
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n n n „ 

9k!,i^ ^ fl*l,fc2,i "S 



+ ■ ■ ■ + 



dgk u -,k\-i 
t) 



dfx J ° d} 2 

we observe that the following induction relations hold: 

(i. /S) 



+ 9k 1 ,...,k x _ 1 ,i 
k 1 ,...,k x _ 1 = l 



+ g, ( h 4r + h 4r + ■ ■ ■ + A 



+ 



9/a 



^i.x (^*ij*ai—)»A-i) 



To obtain the explicit version of the Faa di Bruno in the case of several variables (x 1 , . . . , x n ) 
and one variable y, it suffices to extract from the expression of Y 4l if . provided by Theo- 
rem 3.73 only the terms corresponding to /UiAi + • • • + Pd^d = K, dropping all the X terms. 
After some simplifications and after a translation by means of an elementary dictionary, we ob- 
tain a statement. 

Theorem 3.79. For every integer k > 1 and for every choice of indices i\, . . . ,i K in the 
set {1,2, .. . , n}, the K-th partial derivative of the composite function h = h(x l , . . . , x n ) = 
f(g(x l , . . . ,x n )) with respect to the variables x n , . . . , x 1k may be expressed as an explicit 
polynomial depending on the derivatives of f, on the partial derivatives of g and having integer 
coefficients: 



(3.80) 



d K h 



dx 11 ■ ■ ■ dx lK 



EE E E 

d=l l<Ai<---<A ti <K (ii>l,...,ftd>l /LiiAiH hMdA d =K 

d Xl g 



^/jiH hMd f 

(fyMlH V^d 



E 



n 



d> x d"> l<Vi<fii 



... n 



9 



In this formula, the coset g&i^il'—Af-d, d) was defined in equation (3.71); we have made the 
identification: 

(3.81) {1,...,*} = {1:1:1,... ,1:a*i:Ai, ,d:l:l, . . . ,d:p d :X d }; 

and also, for the sake of clarity, we have restituted the complete (not abbreviated) notation for 
the (partial) derivatives of / and of g. 



§4. Several independent variables and one dependent variable 

4.1. Simplified adapted notations. Assume n = 1 and m > 1, let k e N with k > 1 and 
simply denote the jet variables by (instead of (1.2)): 

(4-2) (x,yj,y{,yi,...,yi) G J£ m . 
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Instead of (1.30), denote the K-th prolongation of a vector field by: 



(4.3) 



£ Y i ^ + £ Y ' 



j=i 



m „ 

y Yi—. 



The induction formulas are: 



(4.4) 



f Y{ :=D 1 (y^-D 1 (X)y{, 
Y{ :=D" (y^-DI {X) V i 



Y{:=D> (YtJ-i? 1 (X)y{, 



where the total differentiation operators D are denoted by (instead of (1.22)): 



(4.5) 



Ox 



£ 



.'7a 



Applying the definitions in the first two lines of (4.4), we compute, we simplify and we organize 
the results in a harmonious way, using in an essential way the Kronecker symbol. Here, the com- 
putations are more elementary than the computations of Y^ and of Y^ ^ achieved thoroughly 
in the previous Section 3, so that we do not provide a Latex track of the details. Firstly and 
secondly: 



m m 

y{ = yi + J2 [y J V H - K *•] yi 1 + E [-K v l ivi> 

^1 =1 h,Z2=l 
m m 

+ E [~ s i Xvw] y^y'y" + E - K 2X »] ^ + 



(4.6) 



Thirdly: 

m m 

ii=i ii,J 2 =i 
(4.7) +E [^'i„'2 tf «3 - S^^i^^] yl 1 yi 2 y[ 3 + 

m 



* 1 > l 2 ^3 j*4 
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+ 



m 

+ E [ 3 ^<Va-< 3^,-^6^]^^ + 
h,h=l 

m m 

+ £ [-<5f 1 3^ i2H! 3-^ 3 3A- !/!l!/l2 ]yi 1 ^^+ [-Si 3 3X yl2 ]y^y l ^ + 

li,fo,l3=i z x , ; 2 =1 

Fourthly: 

m m 
= l Il,f 3 = l 

m 

+ E [^yi V H y i 2v > 3 -€ Q ^v'^]viyM 3 + 

h,h,h=l 

m 

^1 ,^2 ^3 i^4 — 1 

m m 

+ E K Ww] vNNl-vNi 8 + E [ 6 ^ Vl ^ < 4 ^ 3 ] ^+ 

m 

+ E [l23^, lv , a -^6^ Va -^12A; Vl ]„^+ 

m 

(4.8) + E [6^, l8 , 2B ,, -^12^,^-^12^,,]^^ + 

&1,«2,*3=1 
m 

^1 i^4 — 1 

+ E [3^ lHi2 -^12^ i2 ]^^ + 



+ £ 3^, 23/(3 - 8l 2 12X yllyla ] y^y l iy l i + £ [4^ - < 6^] ^ + 

! 1 ,I 2 ,( 3 =l i 1= l 
m 

+ E [ 4 ^, ls , a 4^-^12^]^^ + 

m 



!l,*2,*3=l 



m 

+ E [-< 4^2-C 6 X y i 1 ]&y£+ 

Il,l2 = l 

m m 

+ E - < 4^] yi 1 + E [-< x v* - 5 L 4 -v] ^ 



Ii=l Ii.la=l 

4.9. Inductive elaboration of the general formula. Now we compare the formula (2.9) for 
Y4 with the above formula (4.8) for Y4. The goal is to find the rules of transformation and of 
development by inspecting several instances, in order to devise how to transform and to develope 
the formula (2.25) to several dependent variables. 

First of all, we have to develope the general monomial (y\ 1 ) Ml • ■ • (y\ d . In every monomial 
present in the expressions of Y{, of of Y3 and of Y4 above, we see that the number a of 
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indices lp appearing in all the sums l =1 is exactly equal to /.ti + • • • + [id- To denote 

these fj,x + ■ ■ ■ + indices lp, we shall use the notation: 



(4.10) ha, . . . , h-.m, ■ ■ ■ , Ida, ■ ■ ■ , U-.f. 



MiH K/Jd 



inspired by Convention 3.33. With such a choice of notation, we may avoid long subscripts 

in the indices lp, like ^ fld . It follows that the development of the general monomial 

(y\i Y 1 ■ ■ ■ {v\ d Y d to several dependent variables yields m w+ '" + ' 1 ' ! possible choices: 



\<vi<m i<v d <tid 



where the indices li-i, . . . , lx :fill , . . . , Id-.i, ■ ■ ■ , ld-.fi d take their values in the set {1,2,..., m}. 
Consequently, the general expression of Y£ must be of the form: 



= >*.+£ E E E 

d=l l<Ai<---<A d <« Ml>l,---,Md>l MlAi+--- + Md A d < K + l 
mm mm 

(4- 12 > E ■■■ E E ■■■ E m 

II yiT II 

l<v 1 <f+ 1 l<v d <ii d 

where the term in brackets [?] is still unknown. To determine it, let us examine a few instances. 
According to (4.8) (fourth line), the term [6 y x 2 y — 4 X x 3~\ y~i of Y4 developes as 

E™=i 6 y 3 x i y H ~ S h 4 X x 3 V l 2 in Y 4- Here > 6 y x 2 y just becomes 6 y 3 x2yh ■ Thus, we suspect 

that the term >i j^y!r^n^j^^ • ^-m^i y ^+-+^ d of me second line of (2.25) 

should simply be developed as 



k(k — 1) • • ■ (k — fiiM — ■ ■ ■ — /j,d^d + 1) 
(Ai!)"i m) ■ ■ ■ (Xa^d 

(A \ T,\ v 

v ' ' QK — Vl^l Md^d + MlH \-fj, d yj 



(dx) K ~ tilXl Vd^dQyh.l . . . Qyh-.m . . . Qy'd.l ■ ■ ■ Qy d -Vd 



This rule is confirmed by inspecting all the other monomials of Y{, of of Yg and of Y4. 

It remains to determine how we must develope the term in X appearing in the last two lines 
of (2.25). To begin with, let us rewrite in advance this term in the slightly different shape, 
emphasizing a factorization: 
(4.14) 

«•••(«— /U1A1 - ■ ■ ■ - /i d A d + 2) 



(Ai!)« (Mil ■ ■ ■ (X d \yd Hd\ 



[(/ilAl H + HdXd)X x *- n \ 1 Hd^d + i yfl + '-' + ^d- 1 ] 
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(4.15) 



Then we examine four instances extracted from the complete expression of Y4: 

t m 

m 

£ [-< 4 ^'3-<6^, l!/!2 ]^ 



E H^-Wa^ 



h ^2 ^3>^4~1 



.^1 „ ,^2 „ .^3, ,U 



. *2 '3 



k ! 1 ,! 2 ,!3=1 



and we compare them to the corresponding terms of Y4: 



(4.16) 



[43^3-6*^] ( yi ) s 
[12^2-18^] yi y 
[-10 fy] (yi) 3 2/ 2 , 
[-10A>] ( Vl fy 3 . 



In the development from (4.16) to (4.15), we see that the four integers just before X, namely 
6 = 6, 18 = 6 + 12, 10 = 6 + 4 and 10 = 4 + 6, are split in a certain manner. Also, a single 
Kronecker symbol Sj is added as a factor. What are the rules? 

In the second splitting 18 = 6 + 12, we see that the relative weight of 6 and of 12 is the 
same as the relative weight of 1 and 2 in the splitting 3 = 1 + 2 issued from the lower indices 
of the corresponding monomial y^y 1 ^- Similarly, in the third splitting 10 = 6 + 4, the relative 
weight of 6 and of 4 is the same as the relative weight of 1 + 1 + 1 and of 2 issued from 
the lower indices of the corresponding monomial y^y^ y^y 1 ^- This rule may be confirmed 
by inspecting all the other monomials of Y2, Y 2 , of Y3, Y| and of Y4, Y4. For a general 
K> 1, the splitting of integers just amounts to decompose the sum appearing inside the brackets 
of (4.14) as fiiXx, ■ ■ ■ , l^d^d- In fact, when we wrote (4.14), we emphasized in advance 
the decomposable factor (fiiXi + • • • + fidXd)- 

Next, we have to determine what is the subscript a in the Kronecker symbol Sj . We claim 
that in the four instances (4. 15), the subscript a is intrinsically related to weight splitting. Indeed, 
recall that in the second line of (4.15), the number 6 of the splitting 18 = 6 + 12 is related to 
the number 1 in the splitting 3 = 1 + 2 of the lower indices of the monomial y 1 ^ y 1 ^ . It follows 
that the index l a must be the index l\ of the monomial Similarly, also in the second line 
of (4.15), the number 12 of the splitting 18 = 6 + 12 being related to the number 2 in the splitting 
3 = 1 + 2 of the lower indices of the monomial y 1 ^ y 1 ^ , it follows that the index l a attached to 
the second X term must be the index I2 of the monomial y 2 2 . 

This rule is still ambiguous. Indeed, let us examine the third line of (4.15). We have the 
splitting 10 = 6 + 4, homologous to the splitting of relative weights 5 = (l + l + l)+2in the 
monomial y^y-ly^y^- Of course, it is clear that we must choose the index Z4 for the Kronecker 



symbol associated to the second X term —AX y 



thus obtaining — Sj 4 4 X y i ly i 2y i, 



However, 



since the monomial y^yf^ 3 has three indices l\, I2 and Z3, there arises a question: what index 
l a must we choose for the Kronecker symbol Sj attached to the first X term 6 X y 3\ the index l\, 
the index I2 or the index Z3? 

The answer is simple: any of the three indices l\, I2 or I3 works. Indeed, since the monomial 
ViViy l i is symmetric with respect to all permutations of the set of three indices I2, h}, we 
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have 
(4.17) 

m m 

E {-€ *Xy»M*] y'M'y'M* = £ = 

m 

= £ Hi 6<w B u]^ 2 ^ 

In fact, we have systematically used such symmetries during the intermediate computations (not 
exposed here) which we achieved manually to obtain the final expressions of Yj, of Y J 2 , of Yg 
and of Y 4 . To fix ideas, we have always choosen the first index. Here, the first index is l-y; in the 
first sum of line 9 of (4.8), the first index l a for the second weight 12 is I2. 

This rule may be confirmed by inspecting all the monomials of Y^, of Y|, of Y4 (and also of 
Yg, which we have computed in a manuscript, but not copied in this Latex file). 

From these considerations, we deduce that for the general formula, the weight decomposition 
is simply fi^Xi, . . . ,Hd^d an d that the Kronecker symbol S J a is intrinsically associated to the 
weights. In conclusion, building on inductive reasonings, we have formulated the following 
statement. 

Theorem 4.18. For one independent variable x, for several dependent variables (y , . . . , y m ) 
and for n > 1, the general expression of the coefficient ~Y J K of the prolongation (4.3) of a vector 
field is: 
(4.19) 



v. y: • £ £ £ £ 

d=l l<Xi<--<X d <K ^1 >1,. .. ,fj,d > 1 /^lAiH yfi d \ d <K+l 

^ ^ ^ ^ k(k — !)■■■ (k - /X1A1 + ■ ■ ■ + fi d X d + 2) 

'1:1=1 '1: ( . 1 =1 l d:l- 1 l d: t i d - L 

QK—fiiX-i AidAi+MiH HMd yj 

(k - /X1A1 - ■ ■ ■ - fj, d Xd + 1) 



(9a;) K ~' J i A i Hd*ddy ll -i ■ ■ ■ dy 1: »*i • • • • • ■ dy d ' ^d 

are-HiAi AM A d+AiH ^^<^X 

- S J h i /iiAi 

(<?:r) K ~'' lAl fd^d+igylia . . . Qy l i^ n . . . Qyi d a ■ ■ ■ Qy ld '»d 

QK-fl 1 \ 1 Md^ti+MlH 1- Aid ^ 



- J? MdA d 



(i9x)' c_ ' J i Al vd^d+iQyh-.i . . . Qy l i ■ . ■ dy ld 1 ■ ■ ■ dy ld '< t i 



n vz? 1 n v'z 

l<i/l<Aii l<v d <^d 



Here, the notation dy 1 means that the partial derivative is dropped. 



4.20. Deduction of a multivariate Faa di Bruno formula. Let m e N with m > 1, let y = 
{y 1 , . . . , y m ) e K m , let / = /(y 1 , . . . , y m ) be a C°°-smooth function from K m to K, let x e K 
and let 5 1 = .g 1 (x), . . . , g m = g m (x) be C°° functions from K to K. The goal is to obtain an 
explicit formula for the derivatives, with respect to x, of the composition h := fog, namely 
h{x) := f (g 1 ^), . . . ,g m (x)). For A e N with A > 1, and for j = 1 m, we shall 

abbreviate the derivative by g^ x and similarly for h\. The partial derivatives g y if gy A w iU 
be abbreviated by fli,...,l\ • 
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Appying the chain rule, we may compute: 

hi = E fh 9i , 

!i=l 

m m 

h2 = E f'i,h ^i 1 ffi 2 + 51 ^ - 
!i,i 2 =l Ii=l 

m mm 

^3= 51 fh.h.h g l i g'i g'i + 3 51 fh^yiQ* + 51 /'lflsS 

i 1 ,! 2 ,l3 = l ii,i 2 = l ii = l 



(4.21) 



^4 = 51 /ii^^^s^s^ffi's^ + e 51 fh.i2,hg l i g[ 2 g l 2 3 + 

m mm 

+ 3 51 fh,h92 9 l 2 +4 51] fh,hg l i gs fhgi, 

!i,i 2 =l !l,! 2 =l h=l 

m m 

fa = E fh,h,is,u,is g'i g l i g l i g l i g l i + 10 51 Shmsm g l i g l i g\ 3 g l 2 + 
+ 15 E fh,h,i 3 g l i g l 2 2 g l 2 3 + 10 51 /»i.b,i 3 si 1 #1 2 #e> 3 + 

^l?^2?^3 — 1 iEl 1^2)^3—1 

??i m m 

+ 10 E /ix.i a 32 1 4 2 + 5 E & 34 2 + E fh As 1 • 

ii,! 2 =l ll,I 2 =l (l=l 

Introducing the derivations 
(4.22) 

m m / m Q \ 



^ 3 : = E ^ A + E ai 1 A + E ^ E Ax* ar + E /^^^ J-, 

h = i d 9i i 1= i 9 32 i 1= i \ h =i dih - i 2 , i3 =i d/,2 ' i3 / 

m m Q m a 

F * := g gh J> + g «1 + . . . + g ff h + 

i^i ^i 1 ^ ^ ,fbi ^t: 

m / m m ^ m ^ 

+ E ^ E / ! i.'2^T+ E fh,h,h-Kf 1" 1" E fh,h,-,t\-Kf 

h=l \h=i al ' 2 h,h=i al ' 2 - h h,-,ix=i ' 

w e observe that the following induction relations hold: 

h 2 = F 2 fa) , 
h 3 = F 3 (h 2 ) , 

(4.23) 1 7 ' 



h x = F X (h X -l). 

To obtain the explicit version of the Faa di Bruno in the case of one variable x and several 
variables (y 1 , . . . , y m ), it suffices to extract from the expression of Y J K provided by Theorem 4. 18 
only the terms corresponding to /iiAi + ■ • ■ + PdXd = k, dropping all the X terms. After some 
simplifications and after a translation by means of an elementary dictionary, we may formulate a 
statement. 

Theorem 4.24. For every integer n > I, the n-th partial derivative of the composite function 
h = h(x) = / [g 1 (x), . . . , g m [x]) with respect to x may be expressed as an explicit polynomial 
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depending on the partial derivatives of f, on the derivatives of g and having integer coefficients: 
(4.25) 



dx* ~ ^ 1^ 1^ 2^ (Ai!)^ 11 ui! • ■ • (Ad!)^ d 

d=l l<\^< -<\ d <K n- i >l,...,Hd>l + = K y 1 J ri v a ' ^ a 

m m 

E - E 

11 drM 11 



Qyh.l . . . Qyh-vi . . . Qyl d :l . . . Qy l d:v d -1-1 dx^ 1 dx X 

y y y y l<Vl</il l<v d <Hd 



§5. Several independent variables and several dependent variables 

5.1. Expression of Y? , of Yj 4 and of Y\ 4 4 . Applying the induction (1.31) and working 
out the obtained formulas until they take a perfect shape, we obtain firstly: 



m n m n 

(5.2) Yjj = ^ + E E ^ E E K 5 n ^ 



. '2 



il=lfcl— 1 ^li^2— 1 fcl ,fc2 — 1 



m n 



I ^&3 " 



Secondly: 
(5.3) 

Y U = ^ 2 + E E [^yi^+^ l yin yll -si 1 x% xi2 ]y i k \+ 

1-1=1 fc 1= l 

m 77, 

+ E E [tt^-^^^X'i-^^^Xj^ 4 - 

m ?i 

+ E E [-3,<v2^uR»£«& + 

^l7^2i^3~l ^1 ) ^2 1^3 — 1 

771 77, 

+ E E Kvl 2 ^-^^ 1 ^-^^ 1 ^]^^ 

Z i — 1 fcx ,&2=1 
m tt, 

+ E E ^ [-< « ~ % 

Thirdly: 

y j u . . = ( t + y y k* 1 y>«. „ „ + „ + y* t . „ - a? **a i2 J + 

*1>«2>*3 aj*lx i 2a; £ 3 Z ^ Z ^ |_ *1 x^s'Sy'l *2 x^x^y'l *3 a; ! lx l 2y'l *1 a: 1 a; 2 3- 3 J 17 «1 

Z X =1 fci = l 

m 77, 

+ V V 1 \s kl ' k2 v i +s kl ' k2 v J +s kl ' k2 v j - 

Z_^l Z_^i [ Hi 12 -' X - L 3y l ly l 2 1 »3 > »1 ^x z 2y l ly l 2 >2, '3 ^x z ly l ly l 2 

^1,^2 — 1 fc 1; fc 2 — 1 

<5 fel A" fe2 . , - 5 fcl A" fe . 2 . , -5/ <5 fcl A" fc . 2 . , 1 h! 1 y[ 2 + 

h '1 »'2i'3,'l 22 *2 x'lx l 3y l l l 2 »3 rc'l :E l 2 1 J y *2 1 

m n 

i \ " \ " \ •;k 1 ,k 2 ,k 3 _ s-j srk 1 ,k 2 V fc 3 _ 

Z__, / , [ »1, *2, »3 ^ y l ly l 2 y l 3 I3 *1. *2 x i 3y l ly l 2 

'l^2^3 = l fci,/c2,fe3=l 

_M rfcl,fe2 vfe3 _M rfel,fe2 vfe3 1 ? /l /2 /3 , 



^1^2 = 1 fcl i^2 >fc3 — 1 



m n 
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(5.4) 



m n 

+ \~^L ^lV ill *3 ^y l l y '2 y*3 ] 2/fei ^2 ^4 ^4 + 

m n 

[^1- *2 *Hc i 3y^l ^3- H 'Hc i 2y I l ^~ ^2- *3 *Hc i ly Z l ~~ 

— 1 T & 2 ~ 1 
m n 

' [ 11 <2 ' * 3 y !l »' 2 i;L ' * 2 ' * 3 y' 1 ^' 2 *i< *a> »3 w ^a'i y i2 

£ 1 , i 2 — 1 ^1 ~1 

_jd rfe2.fc3 _ rj £<S2>fc3 V fe l _ J[J X k 2 ,k 3 V fe l _ 

"ll°>l,i! I'Sjl! %%,i 3 x^yh "ll °«a. »3 s*l»'a 

v*3.' c i v>fc2 _ rj 1^3. fc i v fe 2 _ rJ r^ki v fe 2 

% *2 'Vsyh % %. »3 'Vas'l % °»2, »3 ^ly'l 

— P A fcl,fc2 _ X k i< k 2 y k 3 _ £3 :ki,k 2 yk 3 1 li la , 

% °U. »a x^yh °h %> *3 'Vav'i % %> »3 ^nj/'ij yk i yk. 2 ,k 3 ' r 

m n 

" r L °l3°<i.»2. »s y'lj/'a % %. <2> »3 y l iyh % «i . »a. % j/'i B 'a 

^1^2^3 = 1 ^1 ,^2 7^3 7^4 = 1 

_M ifes,*^*! V fc 2 _ XJ ife3. fe li fe 4 V*=2 

i 2 , is y'ly'3 U h n, i2, i 3 y'lyla 

_:i rfei,fc3,fc4 vfca 1 ,/i ,/2 is , 

% %, i2, »3 ^ ! i H i3j y^y^y^,^^ 



m n 

I \ ^ \ " _rj rfci,fe 2 ,fe3 V*3 _ JW rfe2,'=4,'si V k 3 _ p S{ki,ki,k 2 V fe 3 1 ,/l ,,h I 

"•" 2w %"<1,12,«3 jh %°il,*J,<8 j'l <2 »2, <3 a 'l J yk 1 ,k 2 yk 3 ,k 4 ~ l ~ 

1 1 i i 2 ~ 1 1 ! ^-2 ! ^"3 j ^-4 = 1 
m n 

l^ = l_ k^,k2,k 3 — l 
m n 

i \ " \ Si' :k 1 ,k 2 ,k 3 y,k 4 tj cfe4,fci,fe 2 -yk 3 ?j c-k3.k4.kx -yk 2 

+ 2^ 2^ [ d h d n, is, is "Vi _ d <2 d n, ia, is ~ d l 2 °n, ^, is 

_M Sjk 2 ,k 3 ,k 4 yki 1 i! ! 2 

°ii %, < 2 , i 3 H ! 2j yk x yk 2 ,k 3 ,ki 



^1,^2 — 1 ^1 1^2 j^3 !^4 — 1 



5.5. Final synthesis. To obtain the general formula for i , we have to achieve the synthe- 

sis between the two formulas (3.74) and (4.19). We start with (3.74) and we modify it until we 



reach the final formula for Y^ 

H 7 



We have to add the m + ■ ■ ■ + m sums £i 1:1 =i • • • £L=i ' ' ' ^Z-^ d =v 

together with various indices l a . About these indices, the only point which is not obvious may 
be analyzed as follows. 

A permutation a E S^££::+j$* d) yields the list: 

^ <r(l : 1 : 1), . . . , <t(1 : 1 : Ai), . . . a(l : m : 1), . . . , <r(l : ^ : \ x ), . . . 

. . . ,a(d:l:l), . . . ,a(l:l:X d ), . . .cr(d:/x d :l), . . . , <r(d: [i d ■ Ad), 

In the end of the sixth line of (3.74), the last term a(d : fid ■ Ad) of the above list appears as the 
subscript of the upper index k^u-^-xA of the term X M^w^ti) . According to the formal rules of 
Theorem 4.19, we have to multiply the partial derivative of X k " i - d: ^d- x d) by a certain Kronecker 
symbol Sj . The question is: what is the subscript a and how to denote if. 

As explained before the statement of Theorem 4.19, the subscript a is obtained as follows. 
The term a(d : fid ■ Xd) is of the form (e : v e : j e ), for some e with 1 < e < d, for some v e with 
1 < v e < fi e and for some j e with 1 < 7 e < A e . The single pure jet variable 
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appears inside the total monomial 
(5.8) 



1<1^1 <Ml l< y d<Md 

placed at the end of the formula for Y? i (see in advance formula (5.13) below; this total 
monomial generalizes to several dependent variables the total monomial appearing in the last 
line of (3.74)). According to the rule explained before the statement of Theorem 4. 18, the index 
l a must be equal to l e:Ve , since l e . Ms is attached to the monomial (5.7). Coming back to the term 
a{d:(Xd ■ Ad), we shall denote this index by 

(5.9) le:v e — '■ ^7r(e:i/ c :7 c ) =: Ka(d:fj, d :\ d )i 

where the symbol n denotes the projection from the set 

(5.10) {1:1:1,..., l:/ii:Ai, ,d:l:l, . . . ,d:fJ, d :X d } 

to the set 



(5.11) 



{l:l,...,l:/xi,...,d:l,...,d:ju d } 



simply defined by ir{e:u e :-f e ) := (e:u e ). 

In conclusion, by means of this formalism, we may write down the complete generalization of 
Theorems 2.24, 3.73 and 4.18 to several independent variables and several dependent variables 

Theorem 5.12. For j = 1 , . . . , m, for every k > 1 and for every choice of n indices i\, . . . , i K 

in the set {1, 2, . . . , n}, the general expression ofYj i is as follows: 

(5.13) 



d=l l<Ai<-..<A d <K ni>l,...,na>l liiXi + —+HdXi<K+l 



m m 



m m 



E ••• E E ••• E 

!l:l=l (l: Ml =l id:!^ 1 l d:^ d =1 



E 



kl:l:l,~;kl;l:\. =1 k 



E 

l:^ 1 :l,---, k l:^ 1 :\ 1 



E 



R d:l:l....,Kd:l:A 



kj-A i. ,=1 k 



E 

d: ^d :1 ''"' kd: ^d :X d 



E 



V- ^ k a {l:l:l),---, k <r{l: n :\ 1 ),---, k a (d:^ d .X d ) 
V(l).---.V( fll A 1 ).---.V( t , 1 A 1 + ---+ fld A d ) 



)■■■-■ i^d- x d) re g''l ^l+^+Vd^d 
- l MiAjH ht»d A d " 

qk— M1A1 MdA d +/j!H hMd -yj 

g a ,V( (41 xi+- -+ Mti A d +i) . . . o x K(k) dyh-i . . . dy' d: ^d 



E 



E 



j*!o-(l:l:l)>-.-> fe < r(l:J.i:Al)>"-> fe o-(d:(» 1 j:A e [-l) 



^(1).---.V( M1 A 1 ).---.V( M1 A 1 +-+M d ^ d -1) 



CT p,(Ml, A l).--,(M d ,A d ) Tgs MlAi + --+M d A (f -l 
ClAiH hMd A d K 



^K-^iAi Md A d+MlH ^Md^ fc <r(<i:f» (i :A < j) 



| g a; ! x( fll A 1 +-- + ^ £i A (i ) . . . dx'rMQyh:! . . . Qy 1 *v (d : M d : A d ) . . . Q^^^d 



n 



yk 1 , ui:1 ,....k 1 . ul . Xl 



n 



i<"d<Md 



In this formula, the coset S^^^r+u^, was defined in equation (3.71); as in Theo- 
rem 3.73, we have made the identification: 

(5.14) {!,...,«} = {1:1 l:^i :Ai, , «f:l:l, . . . , d:/i d :X d }. 
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5.15. Deduction of the most general multivariate Faa di Bruno formula. Let n e N with 

n > 1, let x = (x 1 , . . . , x n ) e K", let m 6 N with to > 1, let ^ = ^ (a; 1 , . . . , x n ), j = 
1, . . . ,m, be C°°-smooth functions from K" to K m , let y = (y 1 , . . . , y m ) G K m and let / = 
fiy 1 , ■ ■ ■ , y m ) be a C°° function from K m to K. The goal is to obtain an explicit formula for the 
partial derivatives of the composition h := fog, namely 



(5.16) h{x\ .. .,*»):=/ (flV, ••■ > ■ ■ ■ ;2 m 0r\ • ■ • ,0) 



For j = 1 , . . . , to, for A £ N with A > 1 and for arbitrary indices it , . . . , i\ = 1, . . . , n, we shall 
abbreviate the partial derivative g^.f 9 g x ; A by g J i± f and similarly for /i^,...,^. For arbitrary 

indices li, . . . ,l\ = 1, . . . , m, the partial derivative gyf gy A w iU be abbreviated by fi 1: ...,i x . 

Appying the chain rule, we may compute: 
(5.17) 

/in 



ft-il ,i 2 



hiiMM= E fh,h.i 3 [all oil all] + E /fi.b [fflj 5i2,i 3 +ffi2ff!?,i 3 +Si3 3ii,« 2 ] + 

! 1 ,i 2 ,! 3 =l ii,! 2 = l 

m 

+ J2 fh [slU.J 



m 

h h,i2,i 3 M = X! fh,hJ3;h [dll 9% 9il dlt] + 

m 

2^ /li.la.ls |_3i 2 5ia fii,i 4 + 5i 3 Si? fli 2 ,i4 + 5ii S; 2 Si 3 , 



'4. 

h,h,h =1 

+9il all 5i 2 ,i 3 + all all g%M + all all all,^] + 

m 

I x -P ^1 ^2 I ^1 ^2 I ^1 ^2 I I 

~r / ^ J i 1,^2 \9i 1 ,%2 9i 3 ,i 4 < 9i l7 i s 9i 2 t *4 ,M i*3 

m 

I x f ^1 ^2 . l\ i 2 I ^1 ^2 . Z ^ ^2 i 

+ J'i. ! 2 [fix 9i 2: i 3 ,i 4 + 9i 2 5ii,i 3 ,i 4 + 5i 3 ffii,»2,«4 + Si 4 ^H.»2>*3j + 

d,i 2 =l 

m 
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Introducing the derivations 



fcl ^1 =1 

n m ,-. n m 



fci=l li=l °9k 1 fe 1 ,fc 2 = l ! 1= 1 °9k l ,k 2 

m q m ^ \ 

+ I] ft 1 X! a7T+ Y fh,h,h -qt— 

Il=l \i2 = l '" 2 i 2 ,! 3 = l ^'2^3 



/t /ft „ ft /ft i~, 

(5 - 18) P ^ .- sr n h 9 , „h a , , 

r i - — Z^ Z^ »ki,< a Ji + Z^ Z^ »*i,*a-i - 1 r 



fe, =1 I, =1 



q ii ^ ^ -«l,«2,» a ll 



+ H H 5fci,fc 2 ,...,fcA-i,i + 

fci,fc2,...,fc A „ 1 = 1 Zi — 1 ^fci).--)fex— l 

m / m q in q 

+ Y 9^ \Y fa' 1 ' JfT + Y fh,h,h a?— 
!i=l \i 2 = l '" 2 i 2 ,i 3 = l "i'a.ia 

H r Z^ Jh,h,h,~;l\ ~ 



'2.'3.---><A 



we observe that the following induction relations hold: 

^ii,i 2 = (^»i) ' 
^ii,i 2 ,...,i A = (^H,i 2 :---,iA-l) ■ 



To obtain the explicit version of the Faa di Bruno in the case of several variables (cc 1 , . . . , x n ) 
and several variables (y 1 , . . . , y rn ), it suffices to extract from the expression of i provided 

by Theorem 5.12 only the terms corresponding to + ■ ■ ■ + fid^d = K > dropping all the X 
terms. After some simplifications and after a translation by means of an elementary dictionary, 
we obtain the fourth and the most general multivariate Faa di Bruno formula. 

Theorem 5.20. For every integer n > 1 and for every choice of indices i\,...,i K in the set 
{1, 2, . . . , n}, the n-th partial derivative of the composite function 



(5.21) h = h(x\ . . . , x n ) = f (.9 V, . . . , x n ), g m (x\ . . . , x n )) 



with respect to the variables x 11 , . . . , x %K may be expressed as an explicit polynomial depending 
on the partial derivatives of f, on the partial derivatives of the g J and having integer coefficients: 
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(5.22) 



d K h 



dx 11 ■ ■ ■ dx l " 



E E E E 

d—l l<Xi< - <X d <K, /.ii>l,...,/.id>l /xiAiH \-fJ-d^d— k 

QfJ>l~\ h^d j 



E 



E 



, , dy ll l ---dv h ^'i---dy ldl ---dy d: ' 1 oi 



n 

> t<wi</ 

... n 



l<U d <fl d 



d Xd g ld : "d 

Q x i*(,d:v d :l) . . . Q x iv(.d:» d :\ d ) 
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